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IX. Atmospheric Absorption of Heat Radiation by Water Vapour. 


By T. G. Cowurne*. 
|Received September 23, 19449. | 


SUMMARY. 


The absorption of heat radiation by water vapour in the atmosphere 
is considered. It is shown that existing experimental results are all 
consistent with Lorentz collision-broadening. A law is derived for the 
variation of absorption with amount of water vapour; this law differs 
somewhat from that proposed. by Elsasser, indicating a slower extinction 
of radiation than the latter. The correction for pressure variations is 
found to be greater than Elsasser’s normally ; the temperature correction 
agrees with his. 


§ 1. Iv an earlier report, the author has discussed the absorption of 
long-wave (heat) radiation by water-vapour in the atmosphere. The 
quantity of prime interest in that report was the fractional absorption A 
inrelatively narrow wavelength ranges which are, none the less, sufficiently 
wide each to contain several absorption lines. The variation of A with 
the quantity w of absorbing vapour, with the pressure p, and with the 
temperature T, was given special attention. In the present paper, 
further experimental results regarding this variation are collected together, 
and are compared with the results of theoretical investigations. 

The variation of A with w is found to follow a law somewhat different 
from that given by Elsasser (1942). It cannot be expressed by any 
simple analytical formula, and is best exhibited graphically. Plotting 
A directly against w is, however, inconvenient, since dA/dw rapidly falls 
off as w increases ; it is preferable to plot A against log,,w. The latter 
method of plotting, in addition to giving a more compact diagram, also 
has the advantage that an increase in the absorption coefficient in a definite 
ratio throughout the range leaves the form of the resulting curve unaltered, 
and simply displaces it bodily to left or right. This facilitates a comparison 
between the absorptions in wavelength ranges where the mean trans- 
parencies are very different, but where the types of absorption, as measured 
by the relative variations in the absorption coefficient throughout the 
ranges, may be expected to be similar. 

The random effects of the distribution of strong absorption lines on 
the absorption curve are found to be smoothed out only if a frequency- — 
range not much smaller than 25 cm.>1 broad is used. At frequencies 
where the mean transparency does not vary greatly, the range can, 
however, be 100 cm.~! or more broad without any serious change in the 
absorption curve. 


* Communicated by the Author, 
SER. 7, VOL. 41, NO, 313.—FEB. 1950 K 
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§2. The following experimental studies have been made of the variation . 
of A with w. 

(i) The short-wave bands p(0-94u), @(1-134) and W(1-38u) have been 
studied by various authors in order to determine spectroscopically the 
total amount of water vapour present in the atmosphere. In fig. 1, 
curve I., are shown Fowle’s (1912) results for ¥’, giving the absorption 
in a range 0-021u broad at 1:47. His results for ranges of similar 
breadths in the bands p and @ give curves very like that for P’. 

The results obtained by other authors for these bands will not be 
considered in detail. Herzing’s (1937) results referred to wavelength 


log w 


Graphs of absorption A against log w. Curve I., Fowle’s observed results for 
Ww’; curves IT., III. and IV., calculated for the range 100-125 em.-1, with 
line-widths 0-4, 0-2 and 0-1 em.-!. The unit of w is, for curve I. 0-316 cm. 
precipitable water; for curve II., 1:04x10-5em.; for cree TIT 
5:2 10-6 em. ; and for curve IV., 2-6 x 10-§ cm. 


ranges not small compared with the band width; Foster and Foskett 
(1945), on the other hand, appear to have used a very narrow range in 
studying the p-band (they do not state definitely its width, but speak 
of a range only a few Angstroms broad). Hand (1940) compared the 
absorption at the centre of the p-band with that at a neighbouring wave- 
length where the absorption, though much smaller, was not negligible. 
(ii) The 6-3 band has been studied by Fowle and Adel.  Fowle (1917) 
made two bolograms (corresponding to w=0-008 and 0-082 cm. of 
precipitable water) with a 60° rock-salt prism ; his remaining bolograms 
were made with much smaller resolution, using a 15° prism. In the latter 
bolograms, all the bands at wavelengths below 5 were completely smeared 
out ; even beyond 5u there was a considerable mixing of wavelengths 
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which (because of the greater intensity of radiation in the shorter 
wavelengths) resulted in an apparent transfer of the absorption properties 
to bigher wavelengths. We therefore use only the results of the 60° 
prism. These are combined with Adel’s (1941) results for the transmission 
of sunlight through the atmosphere between 5- 5u and I4y. Since 
absorption is almost complete between 5-7 and 7, Adel’s results give 
useful information only beyond Ty. Near 7:'75u, Adel’s curves show the 
absorption of water vapour overlaid by that of N,O, but it is possible to 
obtain a reasonable estimate of that part of the absorption which is due 
to water vapour. 


0) 2 3 
log ur ‘i 


Graphs of absorption A against log w. Curves V., VI. and VII., calculated. for 
the range 225-250 cm.—!, with line-widths 0-4, 0-2 and 0-1 cm.-1; curve 
VIII., derived from Strong’s experimental results at 2lu (the actual 
experimental results are shown by dots). The unit of w is, for curve V., 
1-04 x 10-5 cm.; for curve VI., 5-2 x 10-* cm.; for curve VII., 2-6 x 10-§ cm.; 
and. for curve VIII., 10-5 cm. Clustering round curve V. are shown Fowle 
and Adel’s results for 7-25u (©), 7-5u (x), 7:75w (EJ) and 8u (+); the 
unit of w for these. is, for 7-254, 10-2%cm.; for 7-54, 10-1°?> om. ;. for 
7-754, 10-** em. ; and for 8, 10-°8"> cm. 


Plots of A against log w, using Fowle’s and Adel’s results, gave 
curves for 7:25, 7-5, 7:75 and 8. which were almost identical in form. 
Using different units for w at the four wavelengths, the experimental 
points can be made to cluster round a single curve. They are shown 
in fig. 2 clustered round a curve (curve V.) whose theoretical derivation 
is described below (§ 4). 

These experimental results, like those for the ¥’ band, are open to the 
objection that they refer, for large w, to transmission of sunlight through 
the atmosphere, whose pressure and temperature are non-uniloim, 


K2 
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This objection is, however, probably not serious, since most of the 
absorbing vapour is contained in the lowest km. or two of the atmosphere. 

(iii) In the rotation band there are only isolated experimental results, 
save for a single series of data for 2lu obtained by Strong (1941) 
by residual ray methods. Strong has given results for shorter 
wavelengths also; but where data obtained by other methods are 
available, they are to be preferred. The residual ray method is hard 
to use, especially in regions of intense absorption where, owing to the 
length of the path traversed by radiation in the apparatus, it is difficult 
to estimate the intensity of radiation in the absence of absorption. 
Strong’s data for 2lu are shown in fig. 2 (curve VIII.) ; the zero of 
absorption (not determined by him) is taken to correspond to a transmitted 
radiation equal to 1-95 of his units. 


100 


log («1/w5) 


Graphs of absorption A against log (w/w Curv is C : 
age Yo). Curve IX. is Callendar’s ; 
Pee ROLE oe a & ( ie 0) 1X is Callendar’s ; curve X., 
sasser’s ; curve XI., Se maidt’s ; and curve XII., that derived in the 
present paper. The values adopted for log wy) for the four curves are 
0), 0-5, 0-75 and 1, . 


§ 3. Next will be considered various theoretical and semi-empirical 
laws which have been proposed for the variation of A with w. 

(i) Sutherland and Callendar (1943), following a discussion of CO 
absorption by Callendar (1941), suggested that the infra-red apsoRticen 
of most gases might be taken roughly to follow the empirical formula 


A=wilpains ON SOR My a (1) 
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where Ww, is the value of w for which absorption is half complete. This 
is illustrated in fig. 3, curve IX. The formula gives only a moderate 
representation of Fowle’s data for Y’, and those of Fowle and Adel for 
725-8 ; it indicates tdo rapida variation of A with w below A=50 per cent, 
and too slow a variation above A—=80 per cent. None the less, its 
simplicity makes it convenient for use in approximate work. 

(ii) Elsasser (1938) derived a formula for the variation of A with w, 
assuming that the absorption lines are equal and evenly spaced, and that 
the absorption in each case is given by the Lorentz formula 


ky =hyn-Y{(v—v)® +7}. 2 eee 8) 
He found that, when absorption is complete at the line centres, 
ASCO (Vitbin)*); +4 a cis) a ae bce e(O) 


where 6 is the frequency-spacing between the lines. This formula is 
illustrated in fig. 3, curve X. When w is so small that absorption is 
incomplete at the line centres, the error-function law ceases to be 
applicable ; it asserts that Accw? when w is small, whereas in fact 
Accw then. 

Comparison with the experimental results of Fowle and Adel not only 
confirms the inapplicability of the formula (3) for small w; it also shows 
that the error-function law gives too rapid an increase of A with w when 
A>50 per cent. The reason is to be found in the irregularity in strength 
and spacing of the absorption lines, which permits the existence of certain 
wide gaps between the stronger lines. Such gaps being relatively 
transparent, a much greater increase in w is needed to increase the total 
absorption from a low value to a high than would be required if the lines 
were regularly distributed. 

(iii) Schnaidt (1939) determined the variation of A with w numerically, 
making assumptions. similar to Elsasser’s. A complete blackout of the 
line-centres was not assumed, so that Schnaidt’s results for small 
absorptions effectively complete Elsasser’s. Instead of considering the 
lines as overlapping, Schnaidt supposed each to be cut off at a distance 
from the line-centre equal to five times the line-width 2y, the next line 
beginning there. This assumption is, of course, unjustified; but the 
resultant increase in transparency between the lines is much the same as 
would result from occasional wide gaps between the lines. In fact, 
Schnaidt’s absorption curve (fig. 3, curve XI.) is very like Fowle’s curve 
- for Y over the whole range of A. 

§ 4. ‘To test more exactly the effect of the irregularity of the lines, and of 
the incompleteness of absorption at the line-centres, I have determined. 
“numerically the law of variation of A with w for various ranges in the 
rotation band. The Lorentz formula (2) was assumed, and calculations 
were made with a line-width 2y=0-4 cm.—' in the ranges 100-125, 225-250, 
175-200, 200-225, 300-325 and 325-350 cm.~'; also with line-widths 
of 0-8, 0-2 and 0-1 cm.~! for the first two ranges. The strengths and 
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positions of the absorption lines were taken from the work of Randall, 
Dennison, Ginsberg and Weber (1937). The calculations were mainly 
done in 1943; had they been done more recently, smaller line-widths 
might have been used, and the strengths of absorption lines might have 
been varied to accord with the work of Cross, Hainer and King (1944). 
However, the similarity between the results from the different ranges, 
in which the distribution of line-strengths was very different, indicated 
that a moderate variation in line-strengths cannot seriously affect the 
absorption curves. 

In the calculations with 2y=0-4 cm.~!, the total absorption coefficient 
from all neighbouring lines was first found at intervals of 0-2 cm.? 
throughout the range under consideration. The next step was to find 
the fractions of the range for which the logarithm (to base 10) of the 
absorption coefficient lay in each of successive ranges of magnitude 0-05. 
The contribution to the total absorption from these fractions were then 
determined for values of log ,,w at intervals of 0-25, and were summed. 
Similar methods, with a different frequency-spacing, were used for the 
other line-widths. 

The results for 2y=0-4, 0-2 and 0-1 cm.—! are shown in fig. 1 (curves II., 
III., IV.) for the range 100-125 cm.—!, and in fig. 2 (curves V., VI., VII.) 
for 225-250 cm.-!. The second set of curves shows a slower increase 
of A with w than the first ; this is due to the wide gap containing no 
absorption lines between 223-4 and 244-5 cm.—!, which gives appreciable 
absorption only when w is large: Gaps like this, though usually less 
wide, occur frequently as the wing of the band is approached. The 
absorption curves consequently tend to show a decrease in steepness 
as the wave-number increases, though the curve for 225-250 cm.-1! 
actually represents an extreme form of band-wing absorption; the 
curves for 175-200 and 200-225 cm.~! have practically the same shape 
as that for 100-125 cm.~—!, while those for 300-325 and 325-350 em.-! 
have shapes intermediate between this and that for 225-250 cm.-1. 
A reduction in the line-width also leads to a decreased slope of the absorption 
curve. 

Curves L. and LI. of fig. | are very similar in shape. This is satisfactory, 
for it provides a direct experimental check on the theoretical curves. 
Cornell (1937) found mean line-widths 0-57 and 0-4 cm.~! for the p and 
® bands respectively ; thus the line-width in the ¥’ band is also probably 
about 0-4cem.-1. Hence if the distribution of strong and weak lines has 
the same random characteristics in the Y’ and rotation bands, the 
absorption curves for the two should be much the same. Tbe fact that 
this is indeed so is a striking confirmation of the validity of the Lorentz 
formula (2), used in the calculations. 

The experimental results for 7-25—8y are shown in fig. 2 clustering closely 
round curve V. The agreement with this curve supplies a further 
experimental check on theory. The experimental results, like the 
theoretical curve, represent band-wing absorption (7-25 u is about 215 em.~! 
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from the centre of the 6:3 band). The closeness of agreement between 
the two suggests that the line-width in the 6-3u band is not very different 
from 0-4cem.~1, but little weight can be attached to this suggestion ; 
one would obtain a similar absorption curve, assuming a less extreme 
band-wing absorption but a smaller line-width. 

The line-width at the wing of the rotation band was found by Adel (1947) 
to be 0-22-0-24cem.-'. There is no essential reason why the line-width 
_ Should remain the same right down to very low wave-numbers, but 
observations at centimetre wave-lengths actually show a_ line-width 
which is similar in order of magnitude, though a little smaller (Becker 
and Autler 1946). Taking a uniform line-width 2y=0-24 em.-1, the 
numerical calculations just described can be used to estimate the actual 
absorption in various parts of the rotation band, as well as the run of the 
absorption curve. A comparison of the calculated absorptions with such 
scanty experimental data as exist provides a further check on theory. 

In the range 100-125 cm.~-!, Wright and Randall (1933) give trans- 
mission curves for 80 cm. of damp air, from which I estimate an absorption 
of roughly 55 per cent. The humidity is not quoted, but the quantity 
of absorbing water may be estimated as (72) 10cm. ; for this, the 
calculated absorption is 50+6 per cent. In the ranges 300-325 and 
325-350 cm.~1, Hettner (1918) found 2-5 x 10% em. of precipitable water 
to give 52 per cent and 49 per cent absorption respectively. This was 
at 354° K.; I estimate that 2-5 10-3 cm. at this temperature should 
produce the same absorption, in the two ranges, as about 3x 107? and 
3-4 10-3 cm. at room temperature. The theoretical absorptions of these 
amounts of water are 45 per cent and 43 per cent. In the ranges 200-225 
and 175-200 cm.~1, Kiihne (1933) found 77 per cent and 69 per cent 
absorption to be produced by about 3-3x10-%cm. of water. The 
corresponding theoretical absorptions are 84 per cent and 78 per cent. 

In view of the uncertainties, these figures show fair agreement between 
theory and experiment. The theory is admittedly approximate: the 
line-strengths are somewhat uncertain, particularly in the band-wing, 
where centrifugal stretching of the molecules may be important. 
Nevertheless, a comparison between curve VIII. of fig. 2 and the other 
curves of this figure indicates no reason to believe that the shape of the 
absorption curve at the extreme wing of the band is different from that 
at greater wavelengths. 

It appears, then, that a single absorption curve can be used for regions 
near a band-centre, and a second for band-wings ; because of possible 
differencés in line-width, different curves should be used for the 6-3y 
and rotation bands. However, in atmospheric work complication 1s 
avoided, and remarkably little error is involved, if a single absorption 
curve is used at all wavelengths. Such a curve is shown in fig. 3 
(curve XII.); the material from which it is constructed is given in Table I. 
The table is based on the curves of figs. 1 and 2, adopting a line-width. 
rather above 2y=0-25 cm.-}, and a type of absorption approximating 
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more closely to the band-centre absorption of fig. 1 than to the rather 


extreme band-wing absorption of fig. 2. 


Tables are also given for 


line-widths one-half and one-quarter as large, corresponding to pressures 
of one-half and one-quarter of an atmosphere. 


TABLE I. 


Percentage absorption for given logy, (w/w), 
at pressures of 1, } and } atmosphere. 


—0-2 
—0-4 
—0-6 
—0-8 
—1-0 
—1-2 
—1-6 
—2-0 
—2-4 


Absorption Absorption 
logo (w/o) 1 atm, 4 ane + atm. logy (“/o) 3 atm. pine atm. 

42 421 434 0-2 584 58 574 
34. 354 364 0-4 67 66 65 
yl 29 304 0-6 15 734 724 
204 22 23 0-8 824 81 794 
154 174 19 1-0 88 864 85 
114 134 144 1-2 934 92 Ole 

54 74 £ 1-4 97 96 95 

24 4 54 1-6 99 984 98 

I 2 3 1-8 100 100 994 

TaBLeE II. 


Values of wo, for which absorption is half complete, at room temperature 
and at 220° K. (wy in cm. precipitable water). 


logy9 w 

peas Room PRR te 220° 
5:0 -5-2p —1:3 —1-15 
5:2 5-4 —1-6 —1-55 
5-4 —5-6 —1-9 —1-9 
5-6 —5:8 — 2-35 —2-35 
5:8 -6-05 —2°6 —2-65 
6-05-6-2 —2:3 —2-35 
6-2 —6-35 —1-95 —2-0 
6-35-6:-45 —2:5 — 2-55 
6-45-6-7 — 2-65 —2:7 
6-7 -6-9 —2-2 —2-2 
6-9 -7:-2 —1-85 —1-85 
7:2 -7:-4 —1-45 —1-4 
7-4 —7:6 —1-05 —0-95 
7:6 —7:8 —0:6 —0-4 
7:8 -8:0 —0:15 0-1 


logiy w 
Ser Room temperature 220° 
8-l3u 0-6 0-85? 
13-15 0-02 0-3 % 
15-17 —0:5? 0-0 ? 
17-19 —1:0 —0-35 
19-21 —1-4 —0-9 
21-23 —1:75 —13 
23-25 —2-0 —1-6 
25-27 —2-2 —1-85 
27-30 —2-4 —2-15 
30-34 —2:6 —2:45 
34-40 —2:8 229-7 
40-50 —3-0 — 3) 
50-900 —3-1 —3-15 


Table I. gives A for various values of log.) (w/w), Where wy is, as before, 


the value of w for which A is 50 per cent. 


It needs to be supplemented 


by a further table giving w, for different wavelength ranges; this is 
given, for a pressure of one atmosphere, in Table II., columns (2) and (5). 
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The experimental work on which this table is based is, in general, that 
used in constructing Table II. of the earlier report (Cowling 1943) ; 
beyond 344, however, the numerical calculations described above have 
been used, and the residual ray observations of Strong and Woo (1932) have 
been. disregarded. In the range 8-13u, though values of wy are given, 
the absorption may approximate more closely to continuous absorption 
than to the kind of band absorption to which Table I. applies. The 
curves of figs. 1 and 2 enable one to estimate that the values of logo Wo 
at pressures of 3 and } atmos. respectively are 0-25 and 0-51 greater than 
those at atmospheric pressure. 

§ 5. Pressure effects. In the above it has been assumed, in accordance 
with the Lorentz theory of collision-broadening, that a reduction in 
pressure leads to a reduction in line-width in the same ratio. The Lorentz 
theory is usually accepted as valid for the infra-red absorption of gases 
_ (see, ¢. g., Coggeshall and Saier 1947). Nevertheless, some workers have 
doubted its applicability to water vapour absorption, for the following 
reasons. 

The absorption of water vapour, when between about 10 per cent and 
50 per cent complete, approximately obeys the square-root law A ow! 
(see, e. g., Strong 1941). This law can be given a theoretical justification, 
if the Lorentz formula (2) is assumed to be obeyed, and the amount of 
water present is supposed to be enough to black out the centres of 
absorption lines, but insufficient to produce much absorption in their 
overlapping wings. Making these assumptions, it is found that 


oss yhoo (8) 


where 6 is the frequency-range considered, and the summation is over all 
lines whose centres are in the range. Since y is, on the collision- 
broadening theory, proportional to the pressure p, (4) implies that 
Acp!; that is, a reduction in p produces the same decrease in A as a 
reduction in w in the same ratio. However, experiments show that, 
even for values of w such that A ocw!, A varies much more slowly than p}. 
Falckenberg (1938) found that the effect on the total absorption, integrated 
over all wavelengths, of a k-fold decrease in p, is roughly equivalent to 
that of reducing w in the ratio k}; Schnaidt (1939) obtained a similar 
result from v. Bahr’s (1909) experimental results for absorption near. 2-7. 
Though their results referred to a very limited range of variation of, w, 
this Jaw of pressure variation of the absorption has often been taken as 
generally valid, and has been interpreted as implying that the line-width 
is not proportional to the pressure itself, but to its square root. , 
But it is not actually necessary to abandon the Lorentz theory. The 
numerical work described in the last section enables one to see whence 
the apparent contradiction has arisen. It shows that there are no very 
secure theoretical grounds for using the formula (4). No range of values 
of w exists such that the absorption at line-centres is complete, and that 
in the region of overlapping wings weak ; the wings of strong lines begin 
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to encroach on neighbouring lines before the centres of weak lines are 
blacked out. The calculations for the ranges 100-125 and 225-250 cm.“1, 
with line-width 2y=0-2 em.~!, give absorptions such that A ow, roughly, 
when A is about 25 per cent ; but the actual absorption is then never 
greater than % that given by (4), even when in using this equation some 
allowance is made for overlapping by lumping together all lines whose 
centres are less than 0-5 cm. ! apart. 

The numerical work equally shows that, on the Lorentz theory, the 
pressure variation of A is never as fast as proportional to p!. A reduction 
of 2y from 0-2 to 0-1 cm.~! (equivalent to halving the pressure) leads to 
a decrease in A, in the same two frequency-ranges, which is something 
like in the ratio 0-8 when A lies between 10 and 50 per cent ; the decrease 
is less when A lies outside these limits. The effect on A of a k-fold 
decrease in p is, moreover, not equivalent to that of a \/k-fold decrease 
in w; when A is less than 5 per cent a reduction in p produces little 
change in A, whereas if A exceeds 60 per cent a reduction in p has nearly 
the same effect as a proportionate reduction in w. 

$6. Integral absorptions. These remarks refer to the absorption in 
relatively narrow wavelength ranges. Tables I. and II. can also be used 
to determine the integral absorption of radiation from a black body, 
given by 


cs Oe 


where / is the fraction of the black-body radiation in the range in which 

the absorption is A. If the absorbing vapour is at the same temperature 
as the black body, A is, of course, the emissivity of the vapour, expressed 
as a fraction of the total black-body radiation. 

I have determined the values of A numerically for a black body at 
286° K. at pressures of one, a half, and a quarter atmosphere, - The 
values are somewhat uncertain for large values of w, for which absorption 
in the relatively transparent range 8-13. is important. Absorption 
in this range does not appear to vary with w in the way illustrated by 
Table I. Its run at atmospheric pressure can be estimated from Adel’s 
(1941) diagrams; values at other pressures are uncertain since, as 
Elsasser suggests, the continuous part of the absorption here may be due 
to the wings of far distant absorption lines, with an absorption coefficient 
inversely proportional to the pressure. Such uncertainties affect the 
values of A appreciably when w exceeds a few mm. of precipitable water, 
and so the calculations for the smaller pressures have been restricted to 
small quantities of water. For such values of w, the error in A is not 
likely to exceed 1 per cent or so of the black-body radiation, and even 
when A is greater than a few mm. of water the error in the values at 
atmospheric pressure is not likely to be much lar ger. 


The calculations indicate a pressure dependence of A which is appreciable 
even when A is as small as 2 per cent (this corresponds to A being about 
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fee ae in poe more opaque wavelengths.) A four-fold reduction in 
pressure is equivalent to replacing w by 3w when A=2 per cent, by }w 
wee ey Se and by 4w when A=40 per cent. The caleuleved 
race Or w= 1-2 107% and 2-4 x 103 cm. are compared with Falckenberg’s 
(1938) experimental values in Table III. The two sets of values indicate 


TABLE ITI. 


Calculated integral absorptions, 
compared with Falckenberg’s observed values. 


_w=1-2x 10 em. w=2-4~x 10-3 em. 
A (obs.) A (cale.) A (obs.) A (cale.) 
(per cent) (per cent) (per cent) (per cent) 
1 atm. 7:8 10-8 Peter eran 
+ atm. 6-4 8-7 ee a 
4 atm. 4-5 6-9 8-4 9-6 


xraphs of integral absorption A against log w (w measured in cm.) at pressures 
of 1 (curve (a)), $ (curve (b)), and 4 (curve (¢)) atmosphere. 


pressure changes in the absorption in roughly the same ratio. The actual 
values of the experimental absorption do not agree at all closely with 
theory, but this is not surprising. The experimental results refer to 
absorption of radiation from a source at 100°C. ; as Schnaidt (1938) 
pointed out, they need to be corrected by increasing them by about 
30 per cent ; and if they are corrected thus, the values at | atmos. are 
definitely higher than those found by other workers. 
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The calculated values of A at the different pressures are shown in fig. 4. 
Absorption in the ranges 0-5y and 13-l7p is omitted, | since in the 
atmosphere water-vapour absorption in these ranges is dominated by that 
of CO,; this facilitates a comparison with experiments made with 
atmospheric air and corrected by subtracting the estimated CO, absorption. 
For small values of w (less than 10-?> em.) the calculated A at atmospheric 
pressure is about one-sixth less than the emissivity found by Falckenberg 
(1939), and about the same amount greater than that found by Elsasser 
(1941); it agrees fairly closely with that found by Brooks (1941). For 
large values of w (greater than 1 mm.) the calculated A is less than 
Elsasser’s by two or three per cent of the black-body radiation, and a 
similar amount less than of Brooks : it exceeds by a rather smaller amount 
that of Robinson (1947), derived from what he calls observations of 
minimum radiation, but agrees with the mean of Robinson’s observations, 
including those which he regards as subject to extraneous radiation. In 
view of the uncertainties the agreement of the calculated absorptions 
with experiment is satisfactory. 


§ 7. Atmospheric pressures effects. By reason of the non-uniformity of 
pressure the problem of pressure corrections to atmospheric absorption 
is more subtle than the corresponding problem when the whole absorbing 
column is at the same pressure. Consider, for example, absorption of — 
radiation leaving the earth’s surface in a given direction and traversing 
the whole atmosphere.. The pressure correction to the part of the 
absorption which takes place in the upper atmosphere has to be applied 
to an increment to the absorption. At different pressures the separate 
wavelengths of a range contribute slightly different amounts to a total 
absorption A of assigned amount in that range. Ignoring this difference, 
which is probably unimportant save for very small values of A, the pressure 
correction to an increment of A is equivalent to reducing the amount of 
water vapour in the ratio (dA/dw),, : (dA/dw),, where (dA/dw), is the 
value of dA/dw at atmospheric pressure and (dA/dw),, is that at the pressure 
concerned. The value of the ratio depends on A, the absorption in layers 
below that considered. A similar correction has to be applied in con- 
sidering radiation reaching the ground surface from the upper atmosphere. 

Values of the ratio, deduced from the data of figs. | and 2, are given 
_in Table IV. It will be noted that the pressure correction is greater than 
reducing the quantity of water vapour in the same ratio as p?; save for 
very small absorptions the reduction is more nearly in the same ratio 
as p. Table IV. also contains values of (dA/dw), (dA/dw) giving the 
similar reduction in the effective water amount in considering the integral 
absorption ; the values of this ratio vary less with the total absorption 
than do those of the first ratio. 


§8. Temperature effects, A variation in the temperature T affects the 


absorption in two ways. At given pressure, the line-width varies like 
T-7/10. Temperature changes also affect the populations of the different 
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energy-levels ; the high-energy levels, which contribute most to band- 
wing absorption, are but scantily occupied at low temperatures. A 
decrease in temperature therefore leads to an increase in absorption 
near a band centre (for the rotation band at great wavelengths) but makes 
the absorption cut-off at the wings of the band sharper. 

The only genuinely reliable theoretical method of estimating the 
absorption at low temperatures would be by carrying out calculations 
like those described in §4, for each of a series of wavelength ranges. 
A less complicated approximate method has, however, been used. The 
intensities of lines in the rotation band between 75 and 550 cm.~1 were 
first calculated for a variety of temperatures. As the temperature falls 
off, weak lines are found to be reduced in strength more than the strong 
ones ; thus the gaps between strong lines tend to be more accentuated 
and the absorption tends to possess more of the characteristics of band- 
wing absorption. The effect is not large, however, and is offset by the 
increased line-width ; hence it is reasonable to take A to vary with w 
according to Table I., as at ordinary temperatures. Then, to determine 
the absorption at any wavelength, it is necessary only to know wy. 


TABLE IV. 


Values of (dA/dw),/(dA/dw), and (dA/dw),,/(dA/dw) 
. corresponding to p=4 and + atmosphere. 


ela) eG 
dw)y  \dw/a 2 (=) ; (Fa) 


a? 


(per cent) Zatm. atm. ieee, ZSatm. i atm. 
5 0-82 0-6 5 0-7 0-44 — 
10 0-72 0-46 10 0-65 0-39 
15 0-67 0-4 15 0-62 0-36 
20 0-64 0-36 20 0-60 0-34 
30 0-60 0-33 25 0-58 0-33 
40 0-57 0:3 30 0-56 0-31 
50 0-54 0-28 35 0-55 0-30 
60 0-52 0-265 40 0-55 = 
80 0-5 0-25 


To determine w, the formula (4) is used. This formula indicates 
19 per cent absorption when the true absorption is about 154 per cent. 
Hence w, is about ten times the value of w which (4) indicates to give 
19 per cent absorption. Loh 

In using (4) one is, of course, only following Elsasser (1938), who 
calculated the coefficient of w! in this formula for ranges 25 em.~? broad 
between 75 and 550 cm.~!, at temperatures between 200 and 320° K. 
As noted in my earlier report, Elsasser’s calculations were subject to 
systematic errors; my, own results agree with his (as regards order of 
magnitude) between 75 and 375cm.1, but give values of wy only 


122 T. G. Cowling on the Atmospheric 


one-hundredth as large as his in the range 475-550 cm.~!. Nevertheless; 
my results differ little from his in the ratios which qos piv for the values 
of wy at different temperatures for any given wavelength range, and it is 
therefore unnecessary to give my results in detail. They indicate that, 
for wavelengths greater than 50, the absorption at 220° K. is slightly 
greater than that at room temperature; in the range 40-50 the 
absorptions at the two temperatures are roughly the same; while for 
shorter wavelengths the absorption at 200° K. becomes progressively 
the smaller of the two, until near 20 the value of wy at 220° K. is about 
five times that at room temperature. The estimate at 20m cannot be 
regarded as very reliable ; the estimated value of wy at room temperature 
here is at least twice that indicated by experiment. However, the 
calculated ratios of wy at room temperature and at 200° K. have been 
used to construct a table for w, at 220° K.; this is given in Table II.. 
third and sixth columns. The temperature variation at a frequency v 
in the 6-3 band is taken to be roughly that at a frequency |y—yvo| in 
the rotation band, where vy is the frequency of the 6-3 band centre. 


TABLE V. 


Values of A (a (a) at 220° K. and + atm., compared with (b) values uncorrected 
for the change of Pecnation with temperature, and (c) values at 
room temperature and pressure } atm. 


w (cm.) ROS MO a9 Tig): OR eOR At eso§ Akg es Alar 
percentage A (a) 23 4:5 82 13:6 20-65 28-9 37-4 
percentage A (b) 26 5:3 9-7. 16:0 6924-2 ~33-7 43-6 
percentage A (c) 1-7 3-4 64 LEP 17-8 948 ean 


The temperature variation in w, is large only in the relatively trans- 
parent band wings. In the upper atmosphere, where the total amount 
of absorbing vapour is small, its effect on the integral absorption A is 
masked by the increasing importance, at low temperatures, of long 
wavelengths, where water vapour absorbs strongly. Table V. gives 
the calculated integral absorption at 220° K. of black-body radiation 
at this temperature, both when account is taken of the temperature 
variation of wy, and when the values of wy appropriate to room temperature 
are used uncorrected. The calculations refer in each case to a pressure 
of { atmos. ; the corresponding integral absorptions at room temperature 
are also given for comparison. The increase in A due to the increased 
importance of long wavelengths at low temperatures is about twice the 
decrease due to the increased transparency of the band-wings. 

§ 9. Conclusion. In discussing the absorption of long-wave radiation 
in the atmosphere, drastic simplifying approximations are often necessary. 
It is hoped that the data given in this paper will enable the value of such 
approximations to be assessed. For example, the pressure correction 
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which replaces a quantity w of water vapour at pressure p by a quantity 
w/(p/Pq) at atmospheric pressure p, is in most cases an undercorrection ; 
the true correction is often nearer to replacing w by wp/pq- 

However, it must be emphasized that the approximations justified in 
any particular problem must depend on that problem. As noted in the 
last section, the emissivity of a small quantity of water vapour, measured 
as a fraction of the black-body radiation, increases with decreasing 
temperature. But the radiation reaching the ground from great heights 
does not mirror this increased emissivity. The radiation emitted by 
water at great heights is chiefly emitted in very long wavelengths, which 
are subject to great absorption in the lower levels. On the other hand, 
in considering the radiative equilibrium of the stratosphere the increased 
emissivity at low temperatures should be important. 
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SUMMARY. 

The theory of least squares curve-fitting of non-equidistant data is 
developed in terms of orthogonal polynomials and power moments. 
A modification of the Doolittle scheme is proposed which enables prob- 
able error. calculations to be carried out very simply. The modified 
scheme is illustrated by a specific example. 


§ 1. INTRODUCTION. 
In recent years the use of orthogonal polynomials in the fitting of curves 
to equally-spaced observations has been thoroughly investigated by a 
number of authors (Birge and Shea 1927, Aitken 1932, Davis 1935, 
Kerawala 1941, Fisher 1946, Birge 1947). A comprehensive discussion 
of these investigations has been given by Birge (1947). 

When the observations are not made at equally-spaced intervals, the 
special methods of these authors are no longer applicable. In these 
general cases it is necessary to solve the normal equations algebraically. 
The usual method employed is that of Doolittle (1878). In the Doolittle 
method the calculation of the power-seiies coefficients is quite straight- 
forward, but the probable error calculations (Brunt 1917) are somewhat 
involved. 

Probable error calculations are greatly simplified by the use of orthogonal 
polynomials. Fisher (1946) states that the method of orthogonal poly- 
nomials becomes artificial in the general case where the observations 
are not equally spaced. This is true in the sense that the calculation of 
the orthogonal polynomials T;(~) at each point of observation for use in 
the orthogonal moments 22/T;(~) would require an exorbitant amount 
of time. In fact, as Aitken (1933-34) points out, the only practicable 
method is to base the calculations on the power moments Lyx’. 

Although the direct calculation of orthogonal polynomials and 
orthogonal moments is out of the question, it is highly desirable that the 
fitted curve should be expressed in orthogonal as well as in power-series 
form. The advantages of orthogonal representation are :— 

(i) All curves of degree less than p are obtained in the process of 
obtaining the curve of degree p. 

(ii) The probable error calculations for the power-series coefficients 
and for the fitted values are greatly simplified. 


* Communicated by the Author, 
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The purpose of the present paper is to develop the general theory of 
curve-fitting in terms of orthogonal polynomials and power moments. 
By this means the advantages of the orthogonal method of solution are 
obtained without the prohibitive labour of calculating directly the 
orthogona] moments. 

It is found that the arithmetical procedure is very closely’ related to 
that of the Doolittle scheme. The various terms in the Doolittle scheme 
are shown to be identical with terms appearing in the expansion of the 
orthogonal polynomials. The values §,,(=2T?(x)) and B,, (identicai 
with Birge’s R,,) which form the basis of the probable error calculations 
may in fact be calculated from terms which appear in the Doolittle scheme. 

A variation of the Doolittle scheme is proposed which makes use of 

the advantage arising from orthogonal representation. The scheme is 
extended to include the calculation of probable errors both of the power- 
series coefficients and the fitted values. 
_ Aitken (1932-34) in a series of three papers has given a comprehensive 
treatment of the problem of curve-fitting. In the second of these papers 
he develops a general scheme applicable to the case of non-equidistant 
data. It is doubtful whether Aitken’s scheme would be of very great 
practical use except where the values 2x are fairly small. The scheme, 
unlike that of Doolittle, calculates the coefficients 6;,, and so can be 
_extended to include probable error calculations. It is considered briefly 
in §2. 

Tchebycheff originally obtained his polynomials in the form of continued 
fractions. Isserlis (1927) has discussed the problem of curve-fitting by 
means of Tchebycheff polynomials expressed in this form. Though 
interesting from a mathematical viewpoint, his scheme is too involved 
to yield a satisfactory practical procedure. 


§ 2. GENERAL THEORY. 

Let y(a;), 7=1, 2,...m, represent the n observations which are to be 
fitted to a polynomial of degree p in «, the errors in the observed values 
of the variable x being supposed negligible. The polynomial w,(x) to 
fit these observations is to be chosen so that the sum of the squares of the 
residuals y(x,;)—u(%;) should be a minimum. 


(a) Calculation of Tchebycheff Polynomials. 

The polynomial w,(x) will be expanded in terms of polynomials T,(x) 
of degree j=0 to p, which are orthogonal over the particular set of values 
a, at which the observations were made. The polynomials T;() are the 
orthogonal polynomials of Tchebycheff, satisfying the equations 


EL (e)T m=, Pho S.A) 


i=1 
The equations (1) determine the polynomials T(z) except for a single 
arbitrary constant in each polynomial. This constant 1s accounted for 
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by imposing the additional condition that the coefficient of the leading 
term a in the power-series expansion of T,(x) should be unity. 

In order to derive an explicit expression for the polynomial T,(x), it 
is written as a linear function of the polynomials of lower degree. 


T(x) = 99 + ong gs Ty_a(%) +055 oT y_o(@) +... toy + + (2) 


To find the coefficient «;,, (2) is multiplied by T,(~) and summed over the 
values x; Thus 


IT, (e)T (aw) =2LeIT (a) +05, 2 Tz(x)=0 
and 
tyg = —S0T y (0) DTH). 


This can be written as 


oT —Syx/Sex: ees A : . . ; yy ate (3) 
where ; 
i (4 
S,,= 2a, (x) = 2 Tz(z). se 


Sy=20'T, (x) 
= Lah (who pT y1(%)+ « - - +oy0) 
= +k, 07+ Oh, e194, kab atte hi + &p9Sj0+ = e 5 . (5) 


The values of S,, and «,, can then be calculated in succession from the 

sums of the powers §;,,,9 and the previously calculated values §,,,.and 

Cam (m<k). ‘ 
To express T,(%) as a power series 


Ty (e)=a4 +B, 410 1-E... Begs’. ee rs SERB) 


the coefficients B,,, are determined by expanding each term on the right- 
hand side of (2) in powers of # and collecting terms in x*. Thus 


Bye % jn +%;, 4 Bros, r+, hePR+2, Rr + Hoy, 515-1, ks » (7) 


and the coefficients 8,5, x, B42, ~. «+ - can be built wp in successive steps. 


(b) Coefficients of the Polynomial of Best Fit. 
The polynomial w,(%) is written in the form 
Up(@) = ap T p(%) + Op: Ty_1(@)+ ».+ +4;T,(%)+a, .--. (8) 


where the coefficients a; are to be chosen so that the expression 


Dvp(%;) = Z{y(a;) —U, (x; )}? cSiadiahts Vimenbie) 
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is a minimum. Differentiating the right-hand expression with respect 
to a; to obtain the normal equations, 


PLDC ESCA EA VA EA EAI, SE oS Ra Tatiy 
Then 

Ly (x)T;(%)=a,2'T? (2x) 
or 

yay (VL) 2a) eee. eee) oe u(E) 


The power moments M, of the observations are denned in the usual way : 
MiSs ory(D) x eae Wn MGs eee 12) 
The orthogonal moments M; are defined by the equations 
Mj=2'T,(x)y(a). Ge yas 2 bey tye dB) 
The equation for a, may then be written in the form 
a;=M,/S,,. Ae ACI Sie re (14) 
Now T;(a) in (13) can be expanded by using either (2) or (6). We obtain 
respectively 
y= LY (a) Oty, 5-1 Ty_a(@)+ +» « -g0}/855 
(My tay, 5 Mj, tay; 2M) gt... +05 9MQ}/S,; - (15) 
and 
dg= Ly (x){2 +B; ge I+... +Byof/Sy 
={M,+8, >-1Mj_1+8;,;-2My_2+ sitet Ie + BjoMo}/S,;- Mien eh 6) 


Either of the equations (15) and (16) can be used to calculate a;. 
The polynomial of best fit is most commonly required in the form of 
a power series in 2, 


Uiy(2)= Ppp? Ady pj oe tg 6 © (17) 


Substitution of the expressions T,(a%) given by (6) in (8), and collection of 
terms in x, leads to the equations 


Dg =Uy+By 43,9 4st Bisa staat +++ TBpi%p - - - (18) 


An alternative expression is found by substitution of the values 8; ;, 
from equation (7) and grouping of terms in 4, ,. 


Dyyg= Ay +05 41,9 (Oj41 + Bj 42,541 427 2» By 544%) 
+065 4.0,g(Aj+2 +8543, 542% 187 ++: By 442%) + - « 


- . = Lg 8 
0; +05 4.1,70p, 541 +4425 5 2b + Tp np- + 2) 
La ; 
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Either of the equations (18) and (19) can be used to calculate the coefficients 
be 

pg 
If the power-series expansion of w,(w;) is substituted in (10), the normal 


equations become 
Z{y(%j)— (Oppo tbr e+ « » » +Opp%”)}T,(%)=0 


for all j in the range 0 to p. Using (2), the set of normal equations reduces 
to 

L{y(%j)—(bpo tbpi t+ ».» +Opp%”)}xI=0 
or 


bpoSjo p85 41,0 «+ - +OppS;4p,0=My. .. = . (20) 


This is the form in which the normal equations are generally written. 
Substitution of the calculated coefficients 6 in the equations (20) 
(for j=0 and 1, say) provides a check on the arithmetical calculations. 


(c) Practical Methods of Calculation. 


There are two possible ways of carrying out the calculations ; the first 
using the coefficients «, the second the coefficients 8. In each system the 
05% and the Si are calculated. In the first system (the «-system) the 
coefficients « are used directly, while in the second system (the £-system) 
the coefficients 8 are determined from the « and then used in the remainder 
of the calculations. The appropriate equations for a and 6 are (15) and 
(19) in the «-system, (16) and (18) in the 6-system. 

The «-system is essentially the familiar Gauss-Doolittle method of 
solving the normal equations. Table I. shows a typical Doolittle scheme 
(based on that of Goulden (1939)), with the entries in our notation. 

The principal disadvantage of the a-system is that the terms «,,b,, 
in the expression for b,,; are different for polynomials of different degree 
p, and have to be calculated separately for each degree. In the B-system 
the terms f,,a; in the expression for b,; are the same for polynomials of 
different degree. 

The coefficients B,;, are needed for probable error calculations, and so if 
probable errors are to be calculated it would seem best to evaluate the 
B,, and use the B-system. If a number of different sets of observations 
are to be fitted, the values «; being the same for each set, the B,, can be 
calculated once and for all. It is also possible to calculate the a by the 
y-system (as in the Doolittle scheme), then the values Bites and finally 
the b,,; as in the B-system. The £-system is used in the illustrative example 
in §4. 

Aitken’s scheme (1933-34 b) is essentially a variation of the B-system. 
The scheme is based on a method of solution of simultaneous equations 
developed by Aitken (1932). Since the scheme will hardly be understood 
without a knowledge of this method, a detailed account will not be given 
here, A feature of the scheme is that division is postponed throughout 
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the table; for example, the entry corresponding to 8,,—S39+%19Ss0 
is written as SogSs397—S49S29=So9S21 and the number Sy, is entered in a 
divisor column. There are no non-terminating decimals m the table, 
but the later entries become very large. 

In our notation the entries in the top right-hand line of each section 
of Aitken’s scheme are (when divided by the divisor) the coefficients 
Bj,- ‘Thus, although Aitken does not himself carry out probable error 
calculations, the required values 8,, are displayed in his table. As in 
our B-system, all the polynomials of lower degree are obtained in the 
process of calculating a polynomial of given degree. 


TABLE [. 
The Doolittle Scheme. 

1. Soo S10 Soo S30 Mo 
2. Divide by —So, X19 Og X39 —Ag 
3. Soo S30 Sao M, 
4. Oy x (I) L119 L920 105; 30 a19My 
5. Add Sai Sor Sag M, 
6. Divide by —8,, Cre M4 —Oy 
7 Sao S50 M, 
8. a9 x (1) 29529 &20530 a%zoMy 
9. a1 x (5) Oa380,  %a1Sg1 Ha wy 
10. Add Shy Dae M3 
11. Divide by —S8,, go =A 
12. Seo M, 
13. ag9 X (1) Agz9539 a39My 
14. a3, X (5) 31531 %3iM, 
15. d&g X (10) 32530 32M» 
16. Add Ssa M; 
dig? Divide by —Ss3s —As3 

b33= Os 

Os Wg2)33 Ly 

b31= 21032 ggg My 

D39= ©1931 % 29032 3033 Uo 


eee ems I 


(d) Probable Errors. 
The sum of the squares of the residuals is given by the equation 


Sh = ZY (%;) —Bp(%)}” 


=Di{y(x;)—ay—4,T (4;)— aoe On, 1 )(%;)}", 
and so, using (11), 
Srp = Ly?—By Ah —8 1147 — - - —Sip% 
= Ly?—Mja,—M,a,;— ... ale, Mahe (24) 


= Zy? —M9?/So9— M7?/8i11— ora = MSS pp- 
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2v2 can also be obtained by calculating the various u(x;) directly. As a 
check on the arithmetical calculations it may be desirable to compute 
Sv. by both methods. 

The probable error 7, of an observation y(x;) of unit weight is then given 
by the usual formula 


1,=0-6745{2v2/(n—p—l1)}t.  « - - ss (22) 


To find the probable errors of the coefficients of w,(x), these coefficients 
must be expressed as explicit functions of the observations. Thus 


6, =2{T ,(%)/Sprf -Y(%) 
and hence 


0a, 


Plal=2(ger) lve 


SHB .ic)t. ews Ge in ar 


For the coefficients 6,, 


WD yr 
Oy (x;) 


and, expanding b,,,, as a function of the a; using (18), 
[bp l= 2{T 4 (&%)/SpxtBrevs, ele ea% Sis, maces -#P 
=12{1/SzztBe+t, x/Seoa,ntit «-- +Be/Spp} - + - (24) 


since the cross-products 2'T,(%)T,(~) vanish. 

Because the polynomials T;, are orthogonal, the coefficients a, 
behave as independently observed quantities in probable error calculations. 
For example, formula (24) can be written 


yal =2 (52 ) tule 


» /d0b.,\2 
2 pa. ea LIE 
rDyal= 2 (Get) fag. Wee, eRe a PI 


e) Fitted Values. 


If the fitted values are required at various points w (which may or may 
not be members of the set 2;), they may be calculated by means of the 
power series (17). The probable errors of the fitted values at the points 
x are given by the equations 


1?[ ty (%))=1"[Go]+T?(x)r2[a,]+'T2(a)r2[ao]+ . 
=1p{UYSoo t+ Ti()/S11+T3(@)/Soo+ « « .} 
=15{1/Soo+ (®+B10)?/S1a+ (+B o1%+Bop)?/Sao + +. sf. ory (26) 
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(f) Changes of Origin and Scale. 


It is sometimes necessary to change the origin and scale of the variable 
x. ‘Suppose it is desired to convert the polynomial Up, (x) into a polynomial 
U,(€) in the new variable 


ES Eh tt vice ee tee = (27) 
Then 
Up (E) = Ono +b, (mE —h) +b, o(mE—h)2 + warts 
=Cyp + Cp i€ Cae tees Pee nd Ca lee meee CED 
where 
epg 2 (—h)Ab op ange Set Sed) 
Hi) 


To calculate the probable errors it is necessary to express the coefficients 
as explicit functions of the a,. Thus 


ep 
Cyg=M 2 Vy iy, Pest, ab gO tke OO)! 
kaj 
where 
k 
Vag & (—h) 274) Brea, (31) 
q=i 
and so 
By 2 ore by 
rey) ( 2 See) Dig ie dood B iMG SLE eeitaiy C29) 
hi} 


§ 3. TABLES OF SYMBOLS AND FoRMUL, 
(a) Symbols. 
y(x;) observations. 
U(x) least squares polynomial of degree p which fits the y(x;). 
T;(x) Tchebycheff polynomial of degree j. 


Syo das, 
Six LUT (2). 
Ob i coefficient of T;,(%) in orthogonal expansion of T;,(«). 
Byx coefficient of «* in power-series expansion of T,(x). 

j Lay(x), jth power moment of the observations. 
M; XT; (x)y(x), jth orthogonal moment of the observations. 
a, coefficient of T(x) in orthogonal expansion of w,(x). 
by; coefficient of ad in power-series expansion of w, (2). 
Up residual y(a;)—w,,(%;). 
os probable error of an observation of unit weight. 
r[ ay] probable error of a,. 


r[byy] probable error of by ,. 
r[u,(#)] probable error of fitted value u(x). 
(ah)/m, change of variable. 


C coefficient of € in power-series expansion of w,(&). 


‘pi 
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(b) Formule. 


Sio = 20 
Six =8) 440+ Mn, p19; e-1+%, b-29z, kat ++ + H%RoSjo- 
Heh as Shu! kk 
Byr le a 5 k+IPk+1, Re + k-+2Pk+2, eats ao OG 1Pyon B 
M; = Lyx, 
M; =M, 3 +B, 51M J TP ad -» +BioMy 
=M, p%, 3-1 Mj_ + %,5-2Mj_ its » +» +%j>My 
a; =MyS,,. 
Dns =A + By 41, sri TBj 2,74 4at -- + TBpi%p- 
=O +05 41,50, 941% 42,90p grat ++ +%p)O pp 
Sv == Ly®—Soa8—8, a7... —8,,,05. 
=Ly?—Moa)—Mja,— ... —M,@,. 
= 54y?—M2/S,)—M2/8,,— ... —M2/S,.9- 
r =0-6745{2'v2/(n—p—1)}}. 
[bps] =p (1/55 +5 F41,5/S;aa,g4a+ « - + +B 5j/Spp)? 


] 
tip (tt) =Dyg+By,t+by902-F -. . Dyn”. 
r[U,(%)J= pt /Soo+Ti(®)/811 +T3(#)/Soo+ seas +T5(%)/Spp}?- 
(a) 07-8 ya By eat e +Bio 


§ 4. Mopret Forms anv ILLUSTRATIVE EXAMPLE. 

In carrying out polynomial fitting, it is highly desirable to employ 
some form of tabulation for the various quantities to be calculated. 
“Model Forms’ have been developed for the B-system, and these will 
be illustrated in this section by fitting a seelhe to the observations 
listed in Table II. 


TABLE II. 


A rough graph shows that a quadratic or cubic equation would 
probably provide a satisfactory fit for the points. The calculations will 
be carried through for the cubic; in the f-system the quadratic is 
obtained at the same time. 

The first step is to calculate the sums §,9 of the powers a/ at the points 
of observation for j=0 to 6. This is best done by tabulating the powers 
of « in columns. 
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The second step is to calculate the sums M, of the products ay(a) 
at the points of observation. Again this is best done by tabulating the 
products in columns. The sum of the squares of the observations Dy? 
is calculated at the same time. 

To save space the detailed calculations for these two steps are omitted. 


Model Form A. 


j= —Sin/ Six j= M5 /S5i 


+-0-021875 


S 
10 


X10 


x 
kK 20 
X19 10 


— 0-35 

-+-0-021875 

--478-8775 
~-0-0077 


1 
530 


X19 20 


+478-8775 
-—29-929844 


— 179-6086 
+ 10-4754 


ey] 
S30 
O39 
q 
40 
X19 30 


— 179-6086 
+-11-22554 


4 22807-411 
—3-929 


+0-021875 


X20 
Oo4P 10 


~ 29-929844 
+0-00773 


x 
Sia 


478-8698 8. 


Oo 


Sao 


X29 29 


+0:3531925 | aS, 


169-1332 
-(0-3531925 
+22807-411 


—-14332-729 
—-59°737 


Bl 
S31 
O34 

1 
Ss0 


220% 30 
M1535 


+ 22803-482 
—47-61938 


~--26553-43 
+ 5375-66 
+8054-02 


B20 


%30 
%31P19 
&32P 20 


—29-92211 


+ 11-22554 
—1-04167 
—46-66580 


+0-3531925 | S,. 


8414-945 


++. 1-559576 


Cy 
S60 


309 30 
31531 
32530 


+.1296991-1 
— 2016-2 

— 1085888 
-20467-5 


—36-48195 


4197-4100 
+1-0294 


M, 


4198-4394 


+- 1636-283 
— 1408-134 
+-69-724 


4297-873 


+8909-704. - 


—1716-840 
9291-802 
4-2551-908 


+452-970 


De 


— 47-06855- =Bao 


ah 
Moa 


+0-4143911 | Mia, 


Mja. 


-0-0353981 | Mijas 


+-0-00240151 


4-1-559576 


138-4152 
82-2315 
10-5441 


1-0878 


1 
Sgs 


ay? 


2 
ave 


188619 


233-7294 
95-3142 


13-0827 


2-5386 
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Oj p= —Syx|Sxex 


4-2-94.1250 


+-0-009065 
42-950315 


—1-059186 
4-1-891129 


—0-087612 
+.1-803517 
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Model Form A (Cont.) 


+.0-4143911 


- 0-0125023 
| 0-4268934 


~-0:113036 
+0-313857 


aj—=MG/Syy 


+.0-0353981 


.0-0037453 
+.0-0391434 


1-19 10-* 
0-13 


+0-00240151 


5-30 x 10-8 


1-32 x 10-* 


5-30 x 10-6 


2-69 x 10-3 5-38 x 10-4 


(a) Model Form A. 


The sums S,) of the powers / are entered in the first section of Model 
Form A and the values S;,, %j, (=—Sjz/S,%), and B;;, evaluated. 

The power moments M; are entered in the second section, and the 
orthogonal moments Mj and the orthogonal coefficients a; (=M;/S;,) 
evaluated. On the right-hand side the quantity y* is entered, the sums 
of the squares of the residuals Xv? determined, and the probable error 
rz of an observation of unit weight calculated for the cubic. 

Next the power-series coefficients b,,, are calculated. 'The Model Form 
gives the coefficients of the lower degree polynomials as well as those of 
the cubic, | 

Finally, the probable errors r[{bg;| of the cubic coefficients 63; are 
evaluated. 

The coefficients 6,; were checked, using the normal equations 


So0P30 +8 1093 1 +8 20039 ake S30033 7 M,; | 
» 1030 ae S203 i Sgo0s2 ome S033 =M, ’ 


and as a further check the residuals were found by direct calculation of 
the fitted values. To save space, these check calculations are omitted. 
The cubic of best fit is then 


Ug (2) = (1:804-0-098)-+ (0-314 -+.0-028)a + (0-0391 +0-0027)x2 
++ (0-00240-+0-00054)a:, 
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(b) Model Form B. 
Suppose the fitted values u() at the points 


ves 200 15 10) by OES 4010-15 20, 


and the probable errors of these values, are required. 

The quantities T;(w) at the required points are first calculated as in 
the upper section of Model Form B. These values are then used in the 
middle section to calculate the probable errors r[u(a)]. Finally the 
fitted values are calculated in the lower section. 

All probable error calculations are carried out to slide-rule accuracy 
only. 


Model Form B. 
1-56 |— 47-1) —36-5 
Ta(x) are gt" | | Bao | ‘Ta() 


363 —8000 | 62 +942} —36 | —6480 
—3375 | 351 — 2360 
—1000| 1: si ~-409 

ER) RY | 93) +113 
0 » —386 

+125 ¢ 3 —107 

+-1000 +-649 

+3375 

+8000 


189 x 108 


T3/So0 T3/S33 (r[w(a) ]/r)? 


TS 


0-84 15-6 222 238 
0-47 43 29-4 34:3 
0-208 0-518 0-89 1-68 
0-052 0-005 0-068 0-187 
0 0-106 ° 0-007 0-175 
0-052 0-061 0-177 
0-208 : 2-22 3°16 
0-47 4. 47-2 52-4 
0-84 o{ 327 


bo b,x x ‘ Ug(x) 


1-804 6-278 5-64 : eo Gara 
2.2 41-4 
0-17 -£0-30 
+0-91 +010 
++1-80-£0-10 
4-4-65-£0-10 
+-11-25-£0-42 
+ 23-441-7 
+429 14-2 
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§ 5. Discussion or ILLUSTRATIVE EXAMPLE. 


(a) Arithmetical Calculations. 


The calculations in §4 were carried out on a  hand-operated 
Brunsviga calculating machine. Seven-figure accuracy in multiplication 
and division was obtainable without trouble. 

It will be observed that there is a loss of one significant figure in the 
calculations of S33, dz, @3. In the present example the loss in accuracy 
is unimportant. If the scatter of the observations is small, it may be 
necessary to obtain the coefficients to a greater accuracy. This may be 
done by calculating more significant figures for the earlier quantities in 
Model Form A, particularly «19, %29, Boy and 8,,. These can be calculated 
accurately without trouble, since the number of figures in So, Si9, Sao 
is usually small. 

If the scatter of the residuals is small, it is desirable that the terms 
in the calculation of Lv? should be accurate to the same number of figures 
as Ly*. It may be an advantage to use the alternative forms M2/Soo, 
M;?,/S,,, etc. for these terms. If it is not possible to calculate the terms 
with sufficient accuracy, the individual residuals will have to be calculated. 

Unless the person doing the fitting is an experienced computer, it is 
almost essential that each calculation should be checked before proceeding 
with the next. Whenever possible the check should be done in a 
different way. If desired, a check column can be introduced in the 
calculations of the quantities 8, «, as in the Doolittle scheme. However, 
the check column does not provide a safeguard against small errors, since 
one or more of the values which go to make up the check column entry 
usually are considerably larger than the others. 

In the fitting of curves of higher degree than the third, the arithmetic 
becomes extremely tiresome. Fortunately, in many appliéations, a 
quadratic or cubic provides a satisfactory fit. If it is necessary to fit a 
higher degree polynomial, the Model Forms may be extended to inelude 
the extra terms. 


(b) Choice of the Degree of the Polynomial. 


Birge (1947) has proposed a criterion using the values Sv? (as in Model 
Form A) for deciding on the degree of the most suitable polynomial. 
Increasing the degree from j—1 to 7 at first reduces Xv? by a large factor. 
Soon a value p is reached, such that for values of j greater than p the 
reduction in Xv? is small. The polynomial of degree p then provides the 
best fit for the observations. 

In the illustrative example it would appear that the quadratic would 
give almost as good a fit as the cubic. However, if the individual 
residuals are examined, it will be found that the two observations at 
«“=6-1 and x—6-25 contribute 1:10 to the value 1-45 of Sv?. For the 
quadratic the contributions from these observations amount to 1-09 in 
a total of 2:54. In fact, any least squares polynomial of reasonable 
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degree must pass more or less centrally between these points. It would 
be legitimate to subtract the contributions from these observations before 
comparing the values 2v?. The contributions from the remaining 
observations amount to 1:45 and 0-35 for the quadratic and cubic 
respectively. 

The cubic should then provide an appreciably better fit. Comparison 
of the graphs of the quadratic (dotted line) and cubic (solid line) in fig. 1 
shows that this is indeed so. 


Fig. 1 
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1. SUMMARY. 


This paper describes measurements of the energy, angular and frequency 
distributions of protons and alpha-particles in disintegrations of the silver 
and bromine nuclei of nuclear emulsions, at excitation energies comparable 
with and larger than the binding energy (~700 MeV.). 

Only those charged particles with ionization greater than six times 
minimum have been considered (i.e. greater than that of a 50 MeV. 
proton). The star disintegrations were observed in Ilford C2 plates, 
100 » thick, exposed at Jungfraujoch (3457 m.). 


2. RaTIO OF NUMBERS OF ALPHA-PARTICLES TO PROTONS FROM STARS. 


A QUANTITY of interest in the theory of nuclear disintegrations is the ratio 
(a/p) of the numbers of alpha-particles to protons emitted at different 
excitation energies. Fig. 1 shows the result of measuring, for different 
prong number N, the ratio R of numbers of particles, charge >2, to 
singly charged particles. It has already been shown (Harding, Lattimore 
and Perkins 1949) that the great majority of the singly charged particles 
(~85 per cent) consist of protons. In order to obtain the true «/p ratio 
in AgBr nuclei, however, two main corrections to R are required. For 
N <6, the «/p ratio in AgBr is less than R, because of the contribution 
of the light gelatine atoms in the emulsion (Perkins 1949 a). The latter 
have an «/p ratio of about unity. From the known proportion of stars 
of N <6 due to carbon, oxygen and nitrogen (Harding 1949 a), we thus 
obtain («/P) en, for N=3 as indicated in the graph. Secondly, especially 
for high values of N, particles heavier than alpha-particles are frequently 
emitted from stars (Perkins 1949b). The result of subtracting these 


from R is indicated for stars of 7-11 prongs (N=9), 12-16 prongs (N=14) 
and 17-31 prongs (N=20). From the four true values of («/p), 5, thus 
obtained, it is clear that the «/p ratio is not-a monotonic function of 
excitation energy. The initial increase is easily understandable, since as 
the temperature, and therefore kinetic energy of the emitted particles 
increases, the relative height of the potential barrier to alpha-particles 
and protons becomes less important. As N increases further, «/p appears 
to reach a broad maximum in the region N=10. A maximum in the 


* Communicated by Sir George Thomson. 
+ Now at H, H, Wills Physical Laboratory, University of Bristol, 
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ratio at just this excitation energy (~300 MeV) has been predicted by 
Le Couteur (1949) on the evaporation model. It is to be expected that 
the emission probability for alpha-particles will decrease when the nuclear 


temperature exceeds half the binding energy of a nucleon in an alpha- 
particle. | 


Fig. 1. 
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giving true alpha/proton emission ratio. Vertical lines indicate standard 
deviation. 


3. ENERGY DISTRIBUTION OF SINGLY AND DOUBLY CHARGED PARTICLES. 


(i) Protons. 

The method of obtaining the energy distributions of the singly charged 
particles, predominantly protons, and alpha-particles has already been 
described (Harding, Lattimore and Perkins 1949). It was shown in this 
paper, hereafter referred to as HLP, that the energy distribution of the 
protons was consistent with the evaporation model of the nucleus. This 
model assumes that the energy transferred to the nucleus by the incident 
particle is distributed statistically among the nucleons, by successive 
collisions, before disintegration proceeds ; thus any one nucleon has to 
undergo several collisions before it may acquire sufficient energy to escape 
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from the nucleus (we exclude of course those nucleons “ knocked on ~ 
directly by the star-producing particle, and contributing a high-energy 
“ tail’ to the Maxwellian distribution). 

Fig. 2. 
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Knergy distribution of protons and alpha-particles emitted from silver and 
bromine nuclei at different excitation energies. Protons: (a) 150 MeV. ; 
(b) 400 MeV.:; (c) 850 MeV. Full curves show the calculated evaporation 
distributions. Alpha-particles : (d) 150 MeV. ; (e) 600 MeV. ; (f) 1350 MeV. 


The proton energy distributions, as given in HLP, are reproduced 
in fig. 2 (abc). Kach histogram is taken from a group of stars with the 
mean excitation energies indicated. The latter were calculated by 
assuming the emission of that number of neutrons which would result 
in a stable residual nucleus, The curves were calculated assuming the 
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- free-particle model, with corrections. We see that the thermodynamical 
picture of the disintegration holds good for emission of protons (and 
presumably neutrons), even at excitation energies larger than the total 
nuclear binding energy (~700 MeV.). A more detailed analysis, also 
based on the Fermi-gas model, has recently been given by Le Couteur (1949) 
and by Fujimoto*. 

(ii) Alpha-particles. 

The energy distributions of the alpha-particles are also given in 
fig. 2 (def). In the distribution at 150 MeV. excitation energy (fig. 2d) 
which corresponds to stars of small prong number, the contribution 
from the light emulsion nuclei has been subtracted by a method previously 
described (Perkins 1949a). Fig. 2d, together with evidence described 

__ below, shows that the emission of alpha-particles in small stars may be 
treated as an evaporation process, most of the alpha-particles having 
energies greater. than the barrier potential (10-12 MeV.). As the 
excitation energy increases, the distributions extend to higher energies, 
as do those of the protons, but we also obtain an increasing proportion 
of low-energy alpha-particles. In fact, the most probable energy of the 
alpha-particles decreases as the excitation energy increases. The reduction 
in effective barrier may pethaps be explained by a mean increase in 
nuclear radius due to surface oscillations, as suggested by Bagge (1943). 
Even so, it is not possible to understand why, if all the alpha-particles 
are released in a process of evaporation, there is a reduction in barrier 
for them but not for the protons. We must conclude that a considerable 
fraction, at least, of the alpha-particles are released in a different process. 
This will almost certainly be a preliminary process, when the nucleus 
is very “hot ”’, the reduction in barrier large, and thermal equilibrium 
not yet attained. 


4, ZENITH-ANGLE DISTRIBUTION OF PROTONS AND ALPHA-PARTICLES. 


A second difference in behaviour between protons and alpha-particles 
is revealed in their angular distributions with respect to the vertical. 

_ If the stars are produced by primaries travelling downwards—in what 
follows we assume they are travelling vertically downwards—then the 
emitted particlesshould havea downward component ofmomentum. Let us 
assume that the alpha-particles, protons and neutrons leave the nucleus 
in an evaporation process, that is, the parent nucleus has taken up the 
whole momentum transfer from the primary before it disintegrates. 
Since, in producing a large star, the primary is likely to make several 
collisions in traversing the nucleus, the sideways components of momentum 
transfer tend to cancel out (alternatively, if there is only one very violent 
collision, several mesons may be produced and re-absorbed ; all will be 
strongly directed downwards. This is equivalent to several collisions 
by the primary). We assume therefore that the net momentum transfer 
is vertically downwards, and furthermore that it is numerically equal 

* Y. Fujimoto, private communication, 
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to the energy transfer (i. e. the star is produced by a relativistic primary . 
which also has relativistic energy after leaving the nucleus). All these 
assumptions seem reasonably well borne out by the observations of 
Harding (1949b) on the smaller stars, and we want to investigate their 
validity in the larger ones. 

Let « be the downward velocity of the nucleus as a result of the impact, 
@ the average observed velocity of emission of an evaporating particle. 
We define = \/ (2K/M), where E is the average energy of a particle, M the 
mass. Ifw<%, as is indeed always the case, it can be shown that, assuming 
a Maxwellian energy distribution (HLP), the downward excess of particles 
should constitute a fraction 1/(7/2) .u/# of the total, since the angular 
distribution of the particles in the coordinate system of the nucleus 
will be isotropic. w, as calculated from the mass of the nucleus and the 
energy transfer, can therefore be compared with the uw obtained from 
the observed fractional downward excess of particles. _ 

In the accompanying table, it can be seen that, if we take protons, 
neutrons and alpha-particles together, there is reasonably good agreement 
between the calculated and observed values of u (we assume the fractional 

‘downward excess of the neutrons, and their mean emission velocity, 
equal to that of the protons). In the group of smaller stars (7-11 prongs), 
with excitation energies well below the total binding energy, the ratio of 
downward velocities u,/u,, of alpha-particles to protons does not differ 
significantly from unity, in agreement with the hypothesis of evaporation. 
At higher energies, however, the downward excess of alpha-particles 
becomes extremely small and that of the protons correspondingly larger. 


TABLE I. 

Number of tracks in star .. ve 17-31 12-16 7-11 
Estimated excitation energy (MeV.) 

=momentum transfer (MeV./c) 1350 750 420 
Downward velocity of nucleus=w/c 0-018 0-008 0-0045 
Fractional downward excess of— 

protons .. i < -- (98424)/600 (53+23)/513 (21-+16)/279 

alpha-particles .. x -» (6+18)/322 (12+17)/296 (18+13)/168 


Average downward velocity of— 


(alphas-+-protons-+-neutrons) .. 0-013-+0-004 0-008+0-004 0-007 +0-004 

protons w/c Je i". -. 0-025+0-006 0-014-+0-006 0-009-L0-007 

alphas u,/¢ < a -- 0-002-+40-006 0-003-+0-004 0-006--0-004 
Ratio of downward velocities u,/u, 0-08 +0-24 0-21 40-28 0-68 +0-74 


It is not of course possible to say if the process responsible for this 
behaviour, which appears to set in at excitation energies above 400 MeV., 
is also that which produces the reduction in effective barrier described 
above. A possible explanation, put forward by HLP, is as follows. 
Assume that, in a preliminary process, there is a “‘ boiling off” of alpha- 
particles following intense localized ‘‘ heating ’’ in the top of the nucleus, 
Those travelling upwards escape clean away from the nucleus, whilst those 
travelling downwards are stopped, by collisions and lead to a general 
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_ warming up of the entire nucleus, succeeded by evaporation. The down- 
ward excess of alpha-particles which subsequently evaporate together 
with the protons may thus be cancelled out, and the whole nucleus will 
receive an additional downward momentum from the alphas travelling 
upwards in the preliminary process, thus accounting, in some degree 
at least, for the large downward excess of protons. Poor statistics do 
not allow us to investigate such processes in detail. 

A second hypothesis, which is supported by experimental evidence, 
is that heavy fragments are emitted at high excitation energies, and that 
these are usually unstable and disintegrate spontaneously into alpha- 
particles. The latter then have to overcome a potential barrier much 
lower than that of the parent nucleus. Experimental evidence on the 
distribution of such fragments (Perkins 1949 b) is at present not sufficient 
to throw light on any resulting anomaly in the zenith-angle distribution 
of the alpha-particles. : 

It is to be pointed out that the observations above were made on those 
protons of energies below 50 MeV. However, these would include a 
few not produced in the evaporation process, but knocked-on directly 
by the primary particle. Suppose these knock-ons constitute a fraction f 
of allthe protons. Since they should all travel downwards, they contribute 
approximately an additional fractional downward excess f, if f<l. 
f is difficult to estimate, but it can be seen from the histograms of fig. 1 
that it is not likely to exceed 0-1 in the large stars (i. e. less than half 
the total observed excess for protons). This can to a large extent account 
for the increased downward excess in such stars. 


It can be concluded that both the energy and zenith-angle distributions 
of protons, and presumably neutrons, can be explained very satisfactorily 
on the evaporation picture, even at excitation energies considerably in 
excess of the total nuclear binding energy; but that the corresponding 
distributions for alpha-particles, in stars involving more than 400 MeV. 
energy, are inconsistent with such a picture. The elucidation of the 
mechanism of alpha-particle emission awaits further experimental 
evidence. 
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SUMMARY. 


After a review of the theoretical interpretation of the fine structure 
in X-ray absorption edges, new measurements on the K-edges of different 
compounds of iron at low dispersion are described. From these in con- 
junction with our earlier measurements at high dispersion (Coster and 
Kiestra 1948)+ the following conclusions may be drawn with respect to 
these compounds. 


(a). In a region up to about 70 eV. from the edge the behaviour of the 
absorption coefficient is mainly a property of the atom in question. 


(b). Ina region between 70 and 175 eV. the fine structure is determined 
by the immediate surroundings of the atom. 


(c). In the region above 175 eV. the fine structure depends essentially 
on the crystal lattice. 

This is in accordance with the notion that as the energy of the ejected 
electron increases it penetrates further into the lattice. 


§1. THEORETICAL FOUNDATIONS. 


Ir the absorption coefficient of X-rays is considered in its dependence 
upon the frequency of the radiation, it may in general be stated that the 
absorption decreases with increasing frequency. Every time, however, 
that a critical frequency is surpassed, the absorption suddenly increases. 
This is caused by the fact that at these frequencies a new possibility of 
absorption arises, accompanying the removal of, for instance, a K- or 
L-electron from the atom. The corresponding phenomenon is called the 
K- or L-edge. On the short wavelength side of these edges the variation 
in the absorption is not monotonic, but fluctuates over a distance of 
several hundred volts. The difference between successive maxima and 
minima of the absorption coefficient may easily amount to 10-20 per 
cent. Of this phenomenon, which is called the fine structure of X-ray 
absorption edges, a theory based on wave mechanics was given by 
Kronig (I., 1931; II., 1932) for the case where the absorbing element 
forms part of a crystalline solid. According to this theory the material 
waves, representing an electron which by absorption of an X-ray quantum 
is removed from an inner shell of an atom and traverses the crystal lattice 


* Communicated by the Authors. 
+ In the following part we will indicate this article as J, 
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in a certain direction with a certain energy, suffer diffraction at the 
crystal lattice. The Bragg reflections, occurring at certain wavelengths 
and preventing the electron from entering the crystal (forbidden energy 
zones), even when averaged over all directions of ejection, still leave a 
fluctuation in the absorption probability as function of the wavelength 
of the electron or of the corresponding X-ray frequency. 

Already in his first paper Kronig was able to show that several 
peculiarities of the fine structure of X-ray absorption edges could be 
explained on the basis of his theory. It is for instance evident that 
metals, crystallizing in the same lattice type, should have the same type 
of fine structure. It could indeed be shown experimentally (Coster and 
Veldkamp 1932) that for ejected electrons of the same sort (s- or p- etc. 
electrons) and for atomic numbers not too far apart this was the case. 
Furthermore, the distances of the maxima and minima of absorption 
from a suitable zero point corresponding to the average potential in the 
crystal appeared to be inversely proportional to the square of the lattice 
constant, as was predicted by the theory. If the crystal structure is 
changed by using an alloy instead of the pure element, the fine structure 
changes in the expected way (Coster and Klamer 1934, I.), the same being 
true in those cases where the crystal type is altered by heating (Coster 
and Levi 1939). 

The fine structure of X-ray absorption edges is also observed when 
the absorbing element forms part of a polyatomic molecule (Coster and 
Klamer 1934, II.). For this case too a theory was given by Kronig (1932 
IIl.). By absorption of an X-ray quantum an electron of an mner shell 
is removed with an excess of energy from the atom. The corresponding 
material waves are scattered by the other atoms of the molecule, and 
the resultant wave built up from the primary and the scattered waves 
determines the probability of absorption of the different X-ray frequencies. 
At first sight this method of calculation of fine structure might seem in 
principle to be different from that used in the case of crystalline solids ; 
on closer inspection, however, it is clear that also in the last case, instead 
of regarding the Bragg reflection of the material waves at the crystallo- 
graphic planes, we may equally well determine the resultant wave built 
up by the primary wave and the secondary ones scattered by the 
surrounding atoms. In this process the atoms at short distance play a 
larger part than those farther removed. Coster (1935) was able to show 
that metals crystallizing in a cubic close-packed lattice have the same 
type of fine structure as those crystallizing in a hexagonal close-packed 
lattice, as long as attention is paid only to the relative positions of the 
successive maxima and minima of absorption, while some differences 
are noted when also their form is taken into account. This was associated 
with the fact that in both types of crystal lattices the immediate sur- 
roundings of an atom are formed by the same number of atoms, whereas 
differences in the arrangement of surrounding atoms are only observable 
for larger distances. 
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Kronig’s theory treats the electron liberated from the atom as almost 
free, so that the forces exercised by the lattice may be regarded as small 
perturbations. This implies that the rules stated above are valid only at 
some distance from the zero potential. In the case of metal crystals this 
distance appeared to be more than 70 volts from the edge. Moreover, on 
the experimental side, the fine structure, as it is usually observed, arises 
by the absorption of non-polarized X-rays in a polycrystalline substance 
and consequently the fluctuations in the absorption coefficient which 
still remain represent an effect, averaged over all directions of ejection 
in the crystal. We thus understand that the fine structure of absorption 
edges is a rather rough phenomenon, and it may be regarded as a fortunate 
circumstance that in verifying experimentally the above-mentioned rules 
use was made in the present laboratory of X-ray spectrographs with not 
very high dispersion. Indeed, it is not difficult to ascertain that by using 
too large a dispersion we should easily run the risk of overlooking the 
general features mentioned. 

If, however, it is desired to study the phenomena arising within a short 
distance from the main edge, the dispersion formerly used in this 
laboratory is no longer sufficient. A larger dispersion may be obtained 
firstly by the bent crystal method used by Cauchois (1948) and co-workers 
and S..nner (1941). Some American investigators (Barnes 1933, Beeman 
and Friedman 1939, Bearden and Snyder 1941), using the double crystal 
spectrometer, got very interesting results not obtainable with a smaller 
dispersion. They investigated a region within about 50 volts from the main 
edge, about which Kronig’s theory does not give any information. Indeed, 
to explain the fine structure of X-ray absorption edges, in general, one 
has to bear in mind that in the corresponding transition an electron 
from an inner shell is removed to one of the unoccupied levels in the 
lattice. To calculate the corresponding absorption it is necessary to 
determine the density of states of these unoccupied levels. This can be 
done in two different ways :— 


(le) By starting from the wave function of a free electron. 


(2e) By starting in zero approximation from the wave function of a 
free atom. 


As has already been stated, the first procedure has been employed by 
Kronig, who calculated the perturbation caused by a periodic potential 
that serves to approximate the lattice field. It can only be applied if the 
energy of the electron in the lattice is large compared with the variations 
in the potential. 

The second method has been developed by Wigner and Seitz (1933) 
and Slater (1934). It has been applied by Kruttér (1935) to calculate the 
band structure in the case of Cu for the occupied levels with highest 


energy and the unoccupied levels to a distance of 20 volts from the 
Fermi-limit. 


near the FeK-absorption Edges in the Various Oxides of Fe 147 


The K-absorption edge of Cu, taken by Beeman and Friedman (1939) 
is in good agreement with the theory of Krutter. The agreement is even 
better than could be anticipated from the approximations used. 

In recent times it has been shown that by using the photographic 
method with the single crystal tube spectrograph the same dispersion 
and the same resolving power may be obtained as with the double crystal 
spectrometer (Coster and de Lang 1947). As the dispersion and resolving 
power of the spectrograph in question depend upon its dimensions, it is 
necessary to enlarge these. By constructing a spectrograph with a radius 
about five times as large as that formerly used in this laboratory, we 
(Coster and Kiestra 1948) were actually able to obtain results comparable 
in quality to those of the American authors just mentioned. At the 
same time, by applying some new devices we succeeded in compensating 
for the loss of energy caused by the increased dimensions. 

While hitherto attention has been directed principally to metals, 
especially those of the transition periods, it seemed important to carry 
out an analogous investigation for compounds. There the situation is 
more involved due to the presence of several kinds of atoms, the greater 
complication of the lattice structure and the fact that a given sort of 
atom may have different valencies and even different surroundings in 
the same lattice. Measurements at high dispersion have already been 
published by Yoshida (1941, compounds of Fe, dispersion not very great), 
Sanner (1941, compounds of many transition metals), Beeman and 
Bearden (1942, aqueous solution of salts of Ni, Cu and Zn), Coster and 
Kiestra (1948, oxides of Fe and Mn) and Hanson and Beeman (1949, 
compounds of Mn). ; 

In order to bring out also in the case of compounds the typical difference 
in the behaviour of the fine structure close to and far from the edge, 
we present in this publication results on the oxides of Fe taken at low 
dispersion with the purpose of comparing them with our earlier measure- 
ments on the same compounds at high dispersion. 


§2. EXPERIMENTAL METHOD. 


For these experiments the same spectrograph was used as for the 
investigations at high dispersion in I. The tube between the X-ray tube 
and the crystal and the tube containing the photographic plate holder 
were replaced by others, the distance focus-crystal and also crystal- 
photographic plate being reduced from 115 to 32-5cm. in each case. 
In this way the dispersion decreased in the ratio 3-5 while the region 
investigated increased in the same ratio, 7. ¢. to 400 eV. from the absorption 
edge. 
Instead of tungsten molybdenium was now used as target to avoid the 
tungsten Ly emission line at 1675 X.U. Besides the tungsten hot wire 
was replaced by an oxide-coated platinum cathode. Care was taken that 
the NiKa,, lines were absent as much as possible so that also the fine 
structure in the region of 1660 X.U. could be determined. 
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The analysing crystal was a calcite crystal that remained fixed during 
the exposures. In the beginning the crystal had been carefully ground, 
polished and etched. A time of about 20 hours was then necessary for 
one exposure. As the resolving power of the spectrograph is relatively 
small at these dimensions, it became evident that the crystal does not 
need to have a high resolving power either. Therefore during the last 
exposures a crystal was used that had been ground, polished, but not 
etched, giving the same results in about five hours. This reduction in 
the time of exposure was very favourable, the intensity of the scattered 
radiation being diminished by a factor 4, compared with the former 
exposures. 

Agfa Blau Hart plates were used. These are not as sensitive as Ilford 
Zenith plates but they have a finer grain. On the other hand always at 
least two exposures of the same substance were made and only those 
peculiarities were taken for granted which were present on two (or more 
than two) plates. 

From the photographic plates photometer curves were taken with a 
Moll-photometer, using a beam of light of a height of 3 cm., which is the 
height of the photographed spectrum. In this way the grains of the 
plates do not spoil the curves very much. 

In I. all the photometer curves were changed into intensity curves with 
the aid of a blackening scale for each plate. As can be seen from the 
figs. 5 and 6 in I. on the one hand and figs. 7 and 8 on the other, the 
main form of the curves as a whole is practically the same. For that 
reason in the present article only the photometer curves are given. The 
most interesting feature of these curves is the position and form of the 
maxima and minima of absorption. 


§3. RESULTS. 


In order to bring out the influence of the dispersion on the aspect of 
the fine structure we reproduce in figs. 1 and 2 the results for the K-edge 
found for Fe and published already elsewhere*. Furthermore, we show 
in figs. 3 and 4 the analogous photometer curves as found for various 
oxides of Fe, while fig. 5 gives a comparison at low dispersion of FeF,, 
Fe,0¢4 and FeS. 

As appears from figs. 1 and 2, the diminution in the dispersion causes 
a reduction in resolving power. In fact the minima « and £ in fig. 1 are 
no longer separated in fig. 2. In addition the structure-in the edge has 
become entirely invisible. has 

Fig. 1, embracing only a region of about 70 volts, is identical with the 
curve obtained for metallic iron by Beeman and Friedman (1939) with 


* We call attention to the fact that just as in I. we deviate from the usage 
in the earlier literature by taking increasing absorption in the upward direction 
and increasing energy to the right. In addition for Fe the designation of the 
maxima differs from that of Coster and Veldkamp, an additional maximum 
having been discovered at the higher dispersion. 
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the double crystal spectrometer and entirely of the same type as the 
curves obtained for metallic manganese by Coster and Kiestra (1948) and 
by Hanson and Beeman (1949). An interpretation of this curve with 
which we can fully agree, has been given for the first time by Beeman 
and Friedman on the basis of the theoretical considerations of Krutter. 
(1935) for Cu, already mentioned in the introduction. 

Fig. 2, on the other hand, at a distance of more than 70 volts from 
the main edge to a distance of about 400 volts, shows the typical fine 
structure of the body-centred cubic lattice of which the form is in 
excellent agreement with the theory of Kronig. How strongly this fine 
structure is characteristic for this crystal type has moreover been shown 

by Coster and Veldkamp (1932) and Coster and Klamer (1934, 1 
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Comparison of figs. 3 and 4 shows that here, too, at low dispersion 
the maxima « and f are only separated for FeO. The initial absorption, 
clearly apparent at high dispersion as a small maximum to the left of 
the main edge, is practically unresolved at low dispersion. 

The most striking feature of fig. 4 is that the form of the fine structure 
in the region between 70 and 175 eV. is almost identical for all the oxides : 
two minima with one maximum (D) in between. Fig. 5 shows on the 
other hand that for FeS and FeF, in this region the structure is quite 
different. From this we conclude that the fine structure here is mainly 
due to the O-ions surrounding the Fe-ion. An additional argument is 
that the extrema for the ferrous ion, which is larger, lie closer to the 
edge than those for the smaller ferric ion. The fine structure in this case 
may be compared with that of polyatomic molecules in gases. The result is 
also reminiscent. of that obtained by Coster (1935) in the case of cubic 
face-centered and hexagonal lattices already mentioned in §1. 
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Regarding the fine structure for the oxides above 175 eV. we see that 
the extrema lie in different. places for « Fe,0, and y Fe,0s, while y Fe,Os 
and FeO, behave in an identical way. Now it is known that the «- and 
y-modifications have different crystal structures, while y Fe,O, and 
Fe,0, both are spinels with almost the same lattice constant. From this 
one concludes that the fine structure above 175 eV. is principally deter- 
mined by the lattice type of these compounds. For crystals of metals the 
same is true above 70 eV. 

Summarizing, we may say that as the energy of the ejected electron 
is increased the region of the crystal that plays a réle in determining the 
fine structure is extended more and more. 


The writers are greatly indebted to Professor R. Kronig for his advice 
in the completion of this report. 
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ABSTRACT. 


Solutions of the gravitational field equations of relativity are considered, 
so that the metric is non-static and has spherical symmetry, satisfying 
at the same time the additional equations, ©”, ,—0, where the tensor 
Cijz i8 Weyl’s conformal curvature tensor. The additional equations 
always hold good, whatever be the metric, when R,,=0, but when. R,; 40 
it is shown that the only spherically symmetrical solution of physical 
significance is what is known as Schwarzschild’s internal solution. For. 
the most general gravitational metric of spherical symmetry considered 
here, the condition C,,,,—0 reduces to one equation only, which is one 
of those following from C",,, ,=0; and thus the corresponding space-time 
is conformal to flat space-time even when the latter condition alone is 
satisfied. 


1. INTRODUCTION. 
A GRAVITATIONAL field is distinguished from flat space-time by the 
conditions, ; 
(i.) Rysy=0, for flat space-time, i, a) 
(ii.) Ryi.40, for a gravitational field, 


where R,,,, is the Riemann—Chrystoffel curvature tensor. When (ii.) is 
satisfied the distribution of matter and motion follows from the contracted 
tensor R,;. Weyl’s conformal curvature tensor C,,;, connects the con- 
tracted and original tensors by the equation (Eisenhart 1926) : 


: R 
Crge= Brig + (9x Bie—Ire Bij t+Giz By, —GisRre) + 6 (GraIis—IniYix)» + (2) 
It is well-known that 


Crigz—=9, . . . . . . . . . (3) 


is a necessary and sufficient condition for the space-time to be conformal 
to flat space-time. It.can be immediately verified that the contracted 
tensor, C,;, vanishes identically : 


"Cy, =Cy=0. . . . . . . . (4) 


* Communicated by the Authors. 
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Thus the very structures of the three tensors give ten independent 
components for C,,;,, ten for R,; and 20 for Ryjjx- Moreover, the Bianchi 
identity, viz 


x2 


Rin, +R ig, +R yy, 4220; é . . . . . (5) 
leads to 
Cy, h = sRy;, 5 . 5 A . e e A (6) 
where (Hisenhart 1926) 
Ry=Ry, ;—Ry, rt elGigR, 1—9Gak, i)» wake eae 7) 


in the usual notation. It is clear from (6) that in the unoccupied portions 
of a gravitational field where 


R,=0, Sh eee aE) 


y 


we have also 


V7 ate) 
Noten te vee eth. So -(9} 


(9) is also true when the space-time is conformal to flat space-time, 
that is, when (3) holds good, whatever may be the distribution of matter 
and motion. 

Apart from the considerations outlined here little is known of the 
gravitational significance of the tensor C,,,,._ The following is an attempt 
to investigate it further in relation to distributions of spherical symmetry 
for which the most general metric may be taken as 


ds? = — e*dr? —r?(d62-+-sin® 6 dd?)-+e” dé, ee a) 


where A=A (r,t), v=v(r,t). While (3) implies (9) in the general case 
it is found that, for (10), the latter implies the former. 


2. Tue ConpiTion, C,;,=0. 
For the metric (10) the components of R,,; which are not identically 
zero are 
Ry, Ro», Rss, Rus Rig Ry, 
with 
Kegs ivag Gln? Open ge type 2 ee a (LI) 


They condition (3) implies ten equations, many of which are identically 
satisfied, the rest giving only one independent equation, viz., 


ql —e“)+ dre—*(v' — N’)2+-9e-"(4A—Prv +4 d2)—r2e-*( dy"? dy? JN) =0 
(12) 


bf 


= 
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Here, and in what follows, overhead dashes and dots represent: differentia- 
tions with respect to r and t respectively. For the metric (10), the field 
equations of gravitation (Tolman 1934) give 


paeel 1 
8r'T i= en" (- - z) po Po oe Q A . . . . . : . ° (13) 
Be abd ry hae vin y’—N x ~ rw 
(14) 
AS Aral 1 
ante" (~ a =) ape (15) 
4 1A A 
8rT} e’ = —8nT he = Sit eet: Nar 2 cena ee (16) 
Hence (12) may be expressed as ye ge 
3(1—e~4) + 8arX(T2—TI—TH=0. . 2. 4 (19) 
For a perfect fluid described by the metric (10), 
Ti+Ti=p—p, oo Leeann, haps Coe LL 
and 
: Sar 
e“=)]— z pr’, ee Pai Soe oe a ee en are (19) 


in the usual notation. Thus the condition (3) implies (17) and hence 
(19). Schwarzschild’s internal solution for which p is constant satisfies 
(19) and the corresponding space-time is conformal to flat space-time. 


3. THE Conpition, C",, ,=0. 


Let us now examine the relation (9) for the field given by (10). It 
is obvious from (6) and (7) that, if (9) is true, 


Riz, j— By, 1+ 8(GgB, 1— aR, j)=0. Seeman.) 
For all fields satisfying the equations, 
R,=9, - eee a” 6, leks 
(20) is satisfied and hence (9). For fields oc ak by matter, 
Rg —tiggee sal yg eae ees 
and (20) now becomes 
Ty, 2—-3942, =Ty, ;—I9al, ;. pent Roast Se CES) 
Multiplying both the sides of (23) by g” and summing up, we have 
T} j=. HG Slee Oy Sar Ra Erna 


Thus (23) implies the identities (24) which follow from the field equations 
(22). But the additional content of (23) exercises a severe restriction 
on (10) which we now proceed to investigate, 
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Although there are 64 different equations contained in (23) most o 
them are found to be identities—the independent equations being only 
the following four : 


is 0 1 oT 
re T= To. — Testy? ar vey hy ee a eas) 
ae 
rT = oT gab or ARIE mt erie ee (26) 
7) CL uw L 
ar Tas 30” 5 he athe =T rab oy TD yy— de y'T HT yy, (27) 
7) Oly OL 
y 1a— bes ore ae = NT + he" "AT 4 — FAT Ty + (28) 
For a perfect fluid, 
Dv =(P+ PY» PI ps DI ON dT Wine, Nate) (29) 
and 
gv,» =I, Ss hice mal leant. (OU) 
v, being the flow vector. For (10) the field equations give 
Ve=0) Vg==0; ee pr Le py eit A) nee saree ao} 
and from (29) and (30) we have 
T=p—3p. ea i UK ide Came bach g (3) 
The equations (25) to (28) now become 
==(P-Eolaiesaph, es, sar set at ee (33) 
—(p+p)v'v4=4pr, deen ee Ree aur. (34) 
e [ v'+e'4 ie '(p-+e+ +) | =fellort4or(3i—iy], . (35) 
Lpr(\+v')=Apt+ptip'r),  . . . « . (36) 
respectively. From (33) and (34) it follows that 
(p+ p)vi=pt+ p+3p', eat a LE) 
and, on eliminating v, and v, with the help of (29), we get 
3 2 nA—¥ g 
PS | e+ *(»+e) | Bot eae aN tee ee, (38) 


The problem now is to determine 4 and v in (10) subject to the five field 
equations, (13) to (16), and the four equations, (33) to (36). Actually 
(35) and (36) correspond to /=1 and |=4 respectively in (24) and are thus 
implicit in the field equations themselves. In practice it is convenient 
to use (35) instead of (14). From (37) and (15) we get 


2 
Ee het og he OD) 


e7=1— 
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which-is the same as (19). Hence for (10) the conditions (3) and (9) are 
satisfied by the same value of e~*. (13), (33) and (36) are satisfied by 


4 vy’ 8a 
snivtdet ter ==(1- Fert)... () 


The equations (16), (84) and (39) are mutually consistent. 

The solution is ultimately provided by (39), which gives A, and by (35), 
(38) and (40) which give p, p and vy. All these four equations are found 
to be satisfied by Schwarzschild’s internal solution for which p is constant. 
It is interesting to note that p does not appear explicitly in any one of 
the equations. Either p’ >0 or, as is clear from (38), p’< 0, in which case, 


3 
p+ -(p+p) <9, 


or 
ees 
p’<— = (ptp). 


Stability demands that p’ be zero or negative. These limitations on density 
distributions indicate that only Schwarzschild’s solution (p’=0, p=0) is 
of physical interest under the restriction (9) on (10) when R,,40. If 
p’ <0 and hence 
; 3 
p.S (PP), 

we have 

rT 

M’ir< —12n | pr dr <0, 
0 


where 


a 
M=47 | pr? dr. 
J 0 
This case is of no physical interest since M’ <0. 
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ABSTRACT. 

Fluorescence emission spectra of AgCl, AgBr and AgI, believed to be 
characteristic of the pure solids, have been recorded at liquid nitrogen 
and liquid hydrogen temperatures. The spectra are all in the visible and 
lie just to the long wavelength side of the absorption bands of the 
respective halides. No resolution of the spectra was observed down 
to the lowest temperatures; no pure halide specimen emitted any 
phosphorescence or thermal afterglow. The two different structures of 
AgI capable of existence at ordinary and low temperatures (wurtzite and 
zincblende) gave slightly different emission spectra. 

The temperature at which the fluorescence just became visible was 
about 180° K. for AgCl and AglI, but only 77° K. for AgBr. Excitation 
of the fluorescence could be effected by any wavelength which the halide 
was capable of absorbing. 

All three halides emitted an extra band in the red after being exposed 
to light action at ordinary temperature, or slightly decomposed by heat 
treatment, so as to give a dispersion of free’silver in the halide. Large 
quantities of impurity (~1—3 M. per cent) were required in order to alter 
appreciably the characteristic emission of the basic silver halide. 

The low-temperature fluorescence of the pure halides is believed to be 
due to the direct recombination of excited electrons with positive holes 
with the emission of visible radiation, instead of the emission of thermal 
quanta as at higher temperatures. 


1. INTRODUCTION. 

THis work was prompted by Meidinger’s investigation of the low-tempera- 
ture fluorescence of suspensions of silver chloride, bromide, iodide, and 
bromoiodide in gelatin, prepared in the manner of photographic 
emulsions (Meidinger 1939, 1940). It was uncertain whether the pure 
silver halides would fluoresce in the absence of gelatin : Randall observed 
that silver chloride shows a bright blue-green fluorescence at 77° K., 

but thought that some impurity was responsible. We find that all 
three halides fluoresce strongly in the pure state at a sufficiently low 
temperature, that is, they show the so-called ‘‘ fundamental ” or “* lattice ”’ 
émission. 

-* Communication No. 1294H from the Kodak Research Laboratories. 
+ Now at the University of Helsinki. 
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Until ‘quite recently it was thought that pure solids, apart from. 
manganese, rare-earth, and uranyl compounds having screened incomplete 
shells of electrons, did not show fluorescence (Seitz 1939; Mott and Gurney 
1940). Fluorescence in the pure state, however, has been reported for 
NaCl, KCl, NaBr, KBr, (Kudriavzewa 1934) ZnS, ZnO and CdS (Kroger 
1940), and on theoretical grounds one expects pure substances in general 
to fluoresce at low enough temperatures: what calls for explanation is 
not so much the fact of fluorescence as the mechanism by which it is 
quenched at high temperatures. The main interest of the present work 
lies in the finding that AgCk, AgBr and AgI have to be added to the above 
list. All these substances show a low-temperature fluorescence band 
lying on the long-wavelength flank of the absorption band. 

The present paper (Part I.) is concerned with the individual silver 
halides : Part IT. will deal with mixed silver-halide crystals. The results. 
to be described consist mainly of emission spectra, though some observations 
were made of excitation spectra and the temperature-variation of the 
intensity of emission. Except where otherwise stated, all results were. 
obtained at liquid-nitrogen temperature (77° K.). 


2. SPECIMEN PREPARATION AND EXPERIMENTAL T'ECHNIQUE. 
(i) Specimen preparation. 

The specimens were prepared using distilled water and A.R. grade, or 
purer, reagents, e.g. the AgNO, employed in photographic emulsion 
manufacture which is purer than A.R. grade. In the precipitations. 
0-1 M. solutions were employed, as the washing of the precipitate is. 
inefficient with concentrated solutions (Frankenburger, 1923). During 
the preparation and handling, Kodak Wratten Series “ O” safelights. 
were used for darkroom illumination; occasionally preparation was. 
carried out in darkness. No exposure to actinic light was allowed until 
the specimen was at low temperatures. , 

Usually equivalent amounts of silver nitrate and metal halide solutions. 
were mixed, ¢.g. 50 ml. 0-1 M. KCl and 50 ml. 0-1 M. AgNO, (yielding 
“ equivalent crystals’). “* Excess Ag ” specimens were made by pouring 
the halide solution into an AgNO, solution containing 10 M. per cent more. 
than the equivalent quantity of solute, thus ensuring an excess of silver ion 
throughout the precipitation. The converse procedure was employed to. 
make “excess halide” specimens. All precipitates, particularly the 
“ excess ’’ ones, were thoroughly washed and dried before testing. 

To eliminate the possible influence of inclusions of solution in the- 
precipitate certain samples were shaken in water for some hours. In 
the cases of AgCl and AgBr this yielded large grains of about 1 mm. in 
diameter, not single crystals but hard compact aggregates, which,. 
according to Frankenburger, have no_ inclusions. Further special 
preparations of various specimens will be described later, where: 
appropriate. 
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(i) Haperimental technique. 


The specimens were usually loaded into the lower closed end of a 
length of 3-inch glass tubing, which was immersed in a transparent Pyrex 
or quartz Dewar flask containing liquid nitrogen. They were excited 
by the radiation of a 125 w. mercury high-pressure are transmitted by 
a Wratten 18A filter (3130 A—4047 A.). Excitation by the Hg lines 
transmitted by a 10cm. thickness of a compound solution of 41-5 g. 
CoSO,.7H,O0 and 145 ¢. NiSO,.6-7H,O in 1 litre of distilled water 
(2480 A.-3341 A.) (Bowen, 1935) was found to give the same results. 
In this case specimens were mounted as layers on the lower end of a 
glass strip, a thin layer of gelatin being employed as adhesive. Specimens, 
immersed in the quartz Dewar flask of liquid nitrogen, could then be 
irradiated without the exciting radiation having to pass through glass. 

The fluorescent light was focused on to the slit of a Hilger Medium 
quartz spectrograph (with occasional recourse to the Hilger Constant 
Deviation glass instrument), and the spectrum recorded on Kodak P.1200 
or P.1500 plates. A Wratten No. 2 or 2A filter was mounted at the 
condensing lens to eliminate the ultra-violet radiation scattered in the 
direction of the spectrograph. The glass condensing lens eliminated the 
shorter wavelength ultra-violet when the Bowen filter was employed. 
Direct Hg spectra were always recorded with the fluorescence spectrograms, 
for calibration purposes. With the brightest specimens and a slit width 
of 4mm. an adequate exposure time was 1-2 minutes. Precisely the 
same technique was employed for recording the emission of specimens. 
at liquid hydrogen temperature (20° K.). 


Certain spectra were examined with the aid of a photomultiplier cell. 
The emission of the specimen was allowed to pass right through the quartz 
spectrograph, at the farther side of which it was further dispersed by a 
large replica grating and brought into focus; the linear extension of 
the spectrum was increased by a factor of three. The photomultiplier 
cell with a narrow slit in front of its cathode, was moved along a graduated 
scale in the plane of the spectrum, and current readings were taken at 
various regular intervals. With the help of line-spectrum calibration 
of the scale, and the spectral-sensitivity curve of the photomultiplier cell, 
it was possible to obtain a curve of relative intensity against wavelength 
for the fluorescent light. This was much more convenient than correcting 
for spectral sensitivity of the panchromatic plates employed in the 
photographic recordings. 

The excitation of various Hg lines and groups were observed by dis- 
persing radiation from a high-intensity Hg source with the dispersing 
system employed with the photomultiplier. A small quartz Dewar flask 
containing the specimen in a tube immersed in liquid nitrogen was placed 
in the plane of the Hg spectrum and the effects of separate lines observed 
through a grating, which dispersed at right angles to the direction of 
dispersion of the exciting system. 

NZ 
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3. SILVER CHLORIDE. 
(i) Emission spectra. 
All precipitated specimens at 77° K. emitted a broad structureless 
band extending from 4350 A. to 5800 A. (PI. II., fig. 1, a). The colour 
of the fluorescence was blue-green. Measurements with the photo- 


multiplier arrangement showed that it was a single smooth band with a 
peak at 5000 A. (fig. 1). The extent and position of the band agrees with 


Fig. 1. 
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The spectral distribution of fluorescence of AgCl at low temperatures. Intensity 

in arbitrary energy units : 

Solid curve—specimen at 77° K. 

Dashed curve—specimen at 145° K. 
Ordinates of latter curve adjusted to give same intensity at maximum 
. ” hed ? 
so as to show more clearly the shift of the band to shorter wavelengths 

on increase of the specimen temperature. 


Randall’s results (Randall 1939). That it cannot be ascribed to any 
impurity in the precipitating reagents seems clear from the fact that 
specimens precipitated from eleven different metal halides (e.g. NH,, 
Fe, Co, Ni, Mn, Ca) and several different specimens of AgNO, and Ag,SO, 
gave identical spectra, Some examples are shown in PI. IIL., fig. 1, b, to e. 

To eliminate the possibility that structural defects in the precipitate 
were responsible for the low temperature emission, single crystals of 
AgCl (1-2 mm. in size) were grown from ammonia solution and examined. 
Their low-temperature fluorescence was just as intense as that of the 
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precipitates. It was still possible that defects associated with an excess 
of silver were leading to the observed emission, but this was discounted 
by the fact that annealing the large crystals at 400°C. for five hours in 
an atmosphere of chlorine, and then cooling slowly to room temperature 
(which would also tend to eliminate structural defects) gave specimens 
which were even more strongly fluorescent than before the annealing 
treatment. A specimen prepared by the action of pure, dry chlorine 
on Hilger “ Specpure” silver gave an emission spectrum identical with 
those already recorded. 

If a specimen was melted in the open air, either in a platinum or 
porcelain crucible, so that, in the reducing atmosphere of the Bunsen 
flame, slight decomposition took place, with the formation of free silver, 
the emission was always found to include a band in the red (5800 A. to 
6450 A. Pl. III., fig. 1, fandg). This red band could also be induced by 
a short exposure of the specimen to ultra-violet radiation at ordinary 
temperatures, which led to the formation of “ print-out ”’ silver. Careful 
tests showed that the red emission, which was observed with all halide 
specimens after appropriate preparation or treatment, could only be 
excited by radiation which was absorbed by the halide, 7. ¢. there was 
no specific excitation band of the silver impurity. It was induced more 
readily in an “ excess silver ’ specimen. The red emission is obviously 
associated with the silver excess, but whether it is due to small specks 
or atoms of silver acting as impurity centres, or to the formation of 
F-centres, is at present an open question. Melts showing the red 
fluorescence also gave a weak red phosphorescence and a _ thermal 
afterglow occurring in two bursts, the first having a green and the 
second a red colour. 


(ii) Impurity activation. 

The results of these experiments are important in considering the 
validity of the claim that pure AgCl is fluorescent at low temperatures. 
Mn, Mg, Cu and Pb (as chlorides) were incorporated in the AgCl as 
impurities (by melting in halogen atmosphere), in quantities ranging from 
10-4M. per cent to 5M. per cent. It was found that comparatively 
large quantities of impurity (~1-3 M. per cent) were required in order 
to alter the normal AgCl emission appreciably. Such changes as occurred 
usually consisted of an extension of the emission towards the red; the 
emission spectogram of a specimen containing 5 M. per cent Mg is shown 
as an example (PI. ITLI., fig. 1,2). Specimens containing ~1 M. per cent Mn 
showed a fairly strong red emission band to be ascribed to the Mn impurity 
(Pl. III., fig. 1, i); at 5M. per cent Mn the normal AgCl emission had 
practically been suppressed and replaced by a strong red Mn band 
(5920 A—6450 A.; PI. III, fig. 1, j). 

The specimens which contained ~1 M. per cent Mn were particularly 
interesting in that their colour of emission changed from green to red as 
the intensity of excitation was reduced, and showed a change in the same 
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direction as their temperature was raised. This phenomen is analogous 
to that observed in the case of ZnS+Cu, or ZnS+Mn (Gisolf and Kréger, 
1939), and can be explained by the mechanism put forward by Schén (1942) 
and Klasens (1946). It depends essentially on the fact that the energy 
zap between the ground state of the AgCl centres and the 3p band of 
AgCl is less than the corresponding gap for the Mn centres. An analysis 
based on this and certain other simplifying assumptions (for which 
reference must be made to the original sources) leads to expressions 
which are in qualitative agreement with the phenomena observed. 
Perhaps the most interesting fact to be deduced from the above observations 
is that positive holes are mobile in AgCl at low temperatures, at least 
over a distance of a few lattice spacings. 


(iii) The influence of temperature. 


Measurements with a copper-constantan thermocouple have shown that 
the fluorescence of AgCl commences at a temperature of 180° K. From this 
temperature down to 20° K. the emission band (or rather the peak) moves 
over towards longer wavelengths. At 150° K. the emission appears blue 
in colour, at 77° K. blue-green, and green at 20° K. This is well illustrated 
by the three spectrograms shown in PI. III., fig. 2, a—c, recorded at the 
three specimen temperatures noted above. (The emission at 20° K. was 
recorded with the glass spectrograph and thus the spectrogram is on a 
different scale from the other two.) The shift to shorter wavelengths 
with increasing temperature is also illustrated in fig. 1. 


(iv) Hacittation. 


Tests with the arrangement described in § 2, (ii), have shown that 
AgCl is strongly excited by all Hg lines in the ultra-violet region, but 
only weakly by the Hg 4047 A. line. The absorption edge of AgCl is 
at about 4300 A. and it is thus apparent that the fluorescence of AgCl is 
excited by all wavelengths which it absorbs. The emission excited by 
wavelengths in the range 2480 A.-3130 A. was identical in spectral 
distribution with that excited by the longer ultra-violet wavelengths. 


4. SILVER BROMIDE. 
(i) Emission spectra. 


Ordinary precipitates and open-atmosphere melts of AgBr at 77° K. 
showed a faint emission band in the red extending from 5700 A. to 6450 A. 
(PI. IV., fig. 3, a and 6). This red band was strongest in. the slightly 
decomposed melts and the “‘ excess silver” precipitate. Spectroscopic 
tests with Kodak Infra-Red Extra-Rapid plates showed that it did not 
appreciably extend beyond 6450 A. The red emission was suspected 
of being due to a silver impurity in the various specimens, and to test 
this, a specimen free from silver was prepared by melting some AgBr in 
an atmosphere of bromine and then quenching it in distilled water. It 
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was found to be completely non-fluorescent at 77° K. When it had been 
allowed to “ print out” at ordinary temperatures for a few seconds, 
however, the red emission was again observable at low temperatures. 
This is fairly conclusive evidence that the red emission is due to the 
presence of silver. 

The “ excess bromide ” specimen at 77° K. gave an additional faint 
emission band in the green (PI. IV., fig. 3, c) and it was believed that this 
represented the emission of pure AgBr, in analogy to the blue-green 
‘emission of AgCl, but that a much lower temperature was required in 
‘order to observe it at high intensity. This belief was found to be justified 
when observations were made at a specimen temperature of 20° K. 
AgBr gave a very strong green emission in a single band extending from 
4720 A. to 5375 A. with a peak at 5050 A. This is shown in PI. IV., fig. 3d, 
along with the red “ silver ” band which appeared with quite high intensity 
at 20° K. We believe the green emission band to be characteristic of 
pure AgBr, although it must be admitted that there is not such an over- 
whelming weight of evidence in favour of this claim as in the case of AgCl. 


(il) Impurity activation. 


The most striking result in this section, but one which properly belongs 
to Part L1., was the effect of precipitating small quantities (~5 M. per cent) 
-of AgI in AgBr so as to form mixed crystals. Such specimens showed a 
very strong green emission, the spectrum of which was a single band 
extending from 5150 A. to 6000 A. with a peak at 5500 A. (PI.IV., fig. 4, «). 
"The position agrees with that recorded by Meidinger for the emission of 
silver bromo-iodide emulsion layers. We believe, for reasons to be given 
later, that the emission of such specimens is seated in the iodide ion, and 
that the green band is not a “ sensitized ’” bromide band. 

Evidence supporting the latter contention is provided by the fact that 
inclusions of other mono- and divalent halides, ¢. g., CuCl, CuCl,, PbCl,, 
PbBr,, MgCl,, to the extent of 5M. per cent, had only very weak 
-activating effects. In general they gave specimens showing weak emission 
bands in the red and yellow. (The red bands were not due to a silver 
impurity, as the inclusion of impurity by melting was carried out in a 
‘bromine atmosphere.) Some examples of the spectra are shown in 


Pl. IV., fig. 4, b-d. 


5. SttvER IODIDE. 
(i) Emission spectra. 

The essential features of an AgI emission spectrum are well illustrated 
by P. _., fig. 5, a, the specimen in this case being prepared by the action 
_of iodine vapour on spectroscopically pure silver. The main band extends 
from 4400 A. to 4600 A.; the narrow short wavelength band was of 
‘variable position, but in the case illustrated is situated at 4220 A. Ared 
band, in a similar position to the “silver” bands of AgCl and AgBr, 


164 G. C. Farnell, P. C. Burton and R. Hallama on the 


could be induced by exposure of the AgI at ordinary temperatures. 
Almost all precipitates emitted this red band ; it was particularly strong 
in * excess silver’, and aged, precipitates (Pl. V., fig.'5, b and c¢). Some 
specimens also emitted a faint band in the green (PI. V., fig. 5, ¢ and d), 
which was assigned to some impurity, since it was not present in the 
emission of the purest specimen prepared by the action of iodine on silver, 
or in the emission of the majority of precipitates prepared with alternative 
reagents. , 

A puzzling feature of the first recordings of emission spectra of various. 
AgI specimens was the manner in which the positions of the main band 
and narrow band altered, depending on the conditions of precipitation, 
ce. g., the narrow band was emitted at shorter wavelengths in the case of 
‘excess iodide ”’ specimens (PL. V., fig. 5, ¢). It was later realized that 
these variations were almost certainly due to the fact that at ordinary 
temperatures Ag! is precipitated as a mixture of wurtzite-type (hexagonal) 
and zincblende-type (cubic) structures (Wilsey 1923), and that different 
specimens contained different proportions of the two structures. X-ray 
diffraction patterns confirmed this surmise, and so extensive efforts were 
made to prepare Ag! specimens of one single structure. The methods 
employed for obtaining the wurtzite structure were :— 


(a) Precipitation of AgI dissolved in concentrated HI solution by slow 
addition of ethanol (Helmholz 1935). 


(b) Rapid precipitation of Agl by dilution of a concentrated solution 
of AgI in KI (Bloch and Méller 1931). 


(c) Slow cooling of an AgI melt (Bloch and Méller 1931). 


(d) Slow crystallization of AgI from a dry piperidine solution by 
absorption of atmospheric water vapour by the latter. This method 
yielded soft aggregates of external hexagonal form, but which were 
opaque and easily crumbled to a mass of small crystals. 

For obtaining the zincblende structure, the methods used were :— 


(e) Rapid precipitation of Agl by dilution of a concentrated solution 
of Agl in AgNO, (Bloch and Mdller 1931, Kolkmeijer and Hengel 1934). 


(f) Pulverization of any hexagonal specimen, e.g. pieces of slow cooled 
melt of (c) or the soft aggregates of (d), in an agate mortar and pestle. 

Debye-Scherrer powder diagrams showed that in no case was a specimen 
obtained whose structure was completely one or the other type, but some 
approached 90 per cent structural purity. From these it was possible to 
distinguish the features characteristic of the emissions of the two different 
structural types. They are particularly well illustrated by Pl. _., fig. 6, 
a and 6b which show the emission of a piece of slow cooled melt and that 
of a predominantly cubic precipitate prepared by method (e). The narrow 
band is at a shorter wavelength in the case of the hexagonal specimen, 
while the main band is at longer wavelengths. The latter fact is even 
better illustrated by the two low-exposure spectrograms shown in 
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PL V., fig. 6, ¢ and d of the emission of predominantly cubic and hexagonal 
precipitates. A further comparison of predominantly cubic and hexagonal 
pr BOIS ALGE 1 is made j in PI. V., fig. 6,eand/. It should be noted how intense- 
the red “silver ’’ bands are in the emission of the cubic precipitates ; 
their preparation in all cases involved a heavy excess of silver. 

The ¢hange over in emission characteristics when a specimen is 
Ava from hexagonal to cubic by method (f) is illustrated by 
PIE. 7.. The first spectrogram (Pl. V., fig. 7, a) is that of the 
eee of predominantly hexagonal crystals prepared by method (d), 

_the second and third (Pl. V., fig. 7, 6 and c) those of the same crystals 
after being ground for 10 sae 20 init: respectively. The growth of 
the narrow band at longer wavelengths and the shift of the main band 
towards shorter wavelengths are easily visible. X-ray powder diagrams 
confirmed a general change-over from hexagonal to cubic structure. 

A eareful study of all spectrograms and the X-ray powder diagrams of 
the corresponding specimens showed that all the variations in spectral 
emission exhibited were due to the various proportions of the two different 
structural types in the various specimens. For instance, the specimen 
precipitated in excess iodide solution which gave a narrow band at shorter 
wavelengths (Pl. V., fig. 5, ¢) than did other precipitates, was found to 
be predominantly hexagonal in structure. 

A puzzling feature of the low-temperature emission of AgI is the fact 
that the narrow short wavelength emission band falls within the (low- 
temperature) absorption band of the halide. This is illustrated by 
Pl. VL., fig. 8, a—c which contrasts the emissions of predominantly cubic 
and hexagonal specimens (fig. 8, @ and c) with the absorption spectrum 
of a mixed specimen (fig. 8, 6) (tests showed that the wavelength of the 
absorption edge was sensibly independent of the structure of the Agl 
specimen). We believe that the main emission band of AgI, and the emission 
bands of AgCl and AgBr, are due to the direct recombination of the excited 
electron with a positive hole at low temperatures, this process having a 
relatively low probability at higher temperatures compared with the 
probability of radiationless recombination. At the present moment, 
however, we cannot suggest any mechanism which explains the emission 
of the narrow AgI band, or why it is emitted in a spectral region where 
the absorption coefficient is high (~10° cm.~*) (Schell 1911). 

Its intensity, however, does appear to be controlled by the absorption 
of the halide as it was found that the narrow band was much more intense 
from a specimen which had been sedimented into a compact layer, than 
from the same specimen mounted in a tube as a loose, fine powder. ‘This 
is illustrated by PI. VI., fig. 9,a@and 6. The reason why the self absorption 
should be less in the first case is not at all obvious, however. The intensity 
of the narrow band was particularly high in the case of emissions from 
smooth surfaces of melts, and was much reduced when such specimens 
were broken down with little change in structure. It thus appears ee 
the narrow band emission cannot be assigned to surface ions. 
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(ii) The influence of temperature. 


The fluorescence of hexagonal specimens became visible at about 
195° K. and that of cubic specimens at a slightly lower temperature. At 
77° K. hexagonal specimens showed an appreciably higher fluorescence 
efficiency. PI. VI., fig. 10 compares the emissions of a mixed AgI specimen 
at 77° K. (a) and 20° K. (6), identical exposure conditions being employed 
in each recording. It will be seen that the fluorescence is much more 
intense at the lower temperature, and that the spectrum has shifted and 
extended towards shorter wavelengths. The extension towards shorter 
wavelengths is probably connected with the retraction of the absorption 
edge on reduction of the temperature. The narrow bands do not appear 
to be shifted, and have undergone little change in intensity. 


(ill) Haxcitation. 


All Hg lines up to 4047 A. were found to excite specimens of AgI 
strongly, but not lines beyond this wavelength. The absorption edge of 
AgI at low temperatures is at 4250 A. and thus, as in the case of AgCl, 
we see that fluorescence is stimulated by all wavelengths which are 
capable of being absorbed. This result is indeed to be expected if the 
fluorescence observed is characteristic of the pure solid. 


6. PHOTOGRAPHIC EMULSIONS. 


Some attention was paid to the fluorescence of photographic emulsions 
at low temperatures. The emission spectra corresponded to a simple 
‘superposition of the emissions of their silver-halide contents and gelatin. 
The latter, investigated separately, was found to show a blue fluorescence 
at both ordinary and low temperatures ; the spectrum extended over the 
whole of the visible region. It also showed a moderately strong 
phosphorescence at low temperatures, which was reproduced in the 
emulsions. 

No differences in emission were observed between dye-sensitized and 
unsensitized emulsions of otherwise identical composition. Photographic 
desensitizers were also found to have no appreciable influence on the 
fluorescence. It also appeared that emulsion impurities, such as silver 
sulphide specks on the grain surfaces, and silver-gelatin complexes, eté., 
had no important activating or de-activating effects. 

Meidinger has reported that the efficiency of fluorescence of bromo- 
iodide emulsions (96 per cent AgBr, 4 per cent AgI) depends strongly 
on grain size at a given temperature. At 77° K. a coarse-grain emulsion 
gave a strong fluorescence having a quantum efficiency of about 0-9, a 
series of finer-grain emulsions gave progressively lower efficiences, and a 
Lippmann-type emulsion of very fine grain showed no fluorescence at all, 
We have been unable to confirm this result : we have examined a number 
‘of bromo-iodide emulsions of various grain sizes, including a specially 
prepared series having the same constitution and a wide range of grain 
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sizes, and several transparent Lippmann-type emulsions, and all showed 
about the same high intensity of fluorescence. A speculative explanation 
of Meidinger’s results on fine-grain emulsions would be that he employed 
a method of precipitation that enabled the relatively insoluble Agl to 
‘precipitate before the AgBr with little formation of genuine mixed 
crystals. 

The fact that we recorded identical spectra for the emissions of bromo- 
iodide precipitates and bromo-iodide emulsions (except for the weak 
overlying gelatin emission spectrum in the latter case) shows that the 
-green fluorescence of the latter is not due to activation of the bromo-iodide 
by silver sulphide specks, as has been suggested by Vassy and Vassy (1948) 


7. CONCLUSIONS. 


The experimental results which have been described show that the 
‘silver halides are fluorescent in the pure state at low temperatures. 
Besides the evidence provided by the many alternative preparations of 
the halides, the position of the emission band of each halide on the flank 
ofits absorption band suggests that the fluorescence is of the “‘ fundamental ™ 
or “ lattice’ type. We can explain our results in terms of the following 
simple model. If the system in the excited state (reached directly by 
absorption, or after the electron has spent some time in the conduction 
band) receives sufficient thermal energy to raise it; from the minimum 
of the state to a point where it can cross over to a high vibrational level 
of the ground state, it will accomplish the return to the bottom of the 
ground state with the emission of vibrational quanta. At low temperatures 
however, where the probability per unit time of the above process has 
fallen to a low or negligible value, the excited electron will return directly 
to the ground state (recombination with a positive hole) with the emission 
of radiation. It will be shown in Part II. how the results for the variation 
of fluorescence efficiency of silver bromo-iodide fit in with the simple 
explanation given above. 

We expect that careful investigation will show that many pure crystals 
exhibit a low-temperature fluorescence analogous to that of the silver 
halides. We have observed a fluorescence of the cuprous halides which 
is probably of this type. Our interest was centred on the silver halides 
in the hope that information would be obtained of some use in photo- 
graphic latent-image theory. It is known that a developable latent image 
-can be formed in a photographic emulsion by exposure at or below liquid 
nitrogen temperature ; it is supposed that the photoelectrons released 
during the exposure are trapped mainly in shallow traps, and released 
‘on warming up, when also the ionic movements necessary for latent-image 
formation can take place (Berg 1939). It has been found that in all 
-cases there is a fall in sensitivity at low temperatures as compared with 
that at ordinary temperatures. Some of this fall is no doubt due to the 
-wastage of the energy of the exposure as fluorescence. 
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SUMMARY. 

A new method of measuring the scattering of charged particles in their 
passage through photographic emulsions has been developed. It has 
been employed to determine the mass of charged particles emerging from 
nuclear ‘‘ explosions’ produced by cosmic radiation. Of the mesons 
among the “shower” particles, produced in these explosions with 
momentum less than about 180 MeV./c, it is shown that at least 90 per cent 
are 7-particles ; and that less than 5 per cent, and possibly none, are 
p-particles, electrons or other particles of small mass. If mesons of 
greater mass are sometimes created in nuclear interactions with an energy 
less than 300 MeV., they must be produced with a frequency less than 
1 per cent of the protons and less than 5 per cent of the 7z-particles. 
The observations therefore prove that the great majority of the charged 
particles created in nuclear interactions are z-particles. 

The greatly increased rate of measurement with the new method and 
the fac’ that it allows the determination of the sca’tering of particles of 
much greater energy than has been possible hitherto will enable a detailed 
study of the nuclear transmutations produced by protons and «-particles 
with energies in the interval from 10° to 10! eV. to be made. 


INTRODUCTION. 

In Part II. of the present series of papers (Camerini et al. ; 1949), which 
will be referred to as (II.), it was shown that, in favourable cases, it is 
possible to determine the rest-mass and the energy of fast-charged 
particles, emitted during the explosive disintegration of nuclei, by 
observations of the grain-density in the tracks and the deviations in them 
due to multiple Coulomb scattering. The essential principles of the 
methods employed for measuring the scattering have been described by 
Goldschmidt—Clermont et al. (1948) and by Davies et al. (1949). 


* Communicated by Professor ©. F. Powell, F.R.S.+ 

+ It had been intended that Part III. of the present series of papers should 
contain a discussion of the experimental material-given in Parts I. and IT. 
The rapid development of the methods discussed in Part IT.—and, in particular, 
the improvements in the scattering method described by Mr. Fowler in the 
present paper—allowed us, however, to obtain a more detailed insight into 
many features of the nuclear disintegrations produced by protons and a-particles 
of great energy. It therefore seemed desirable to make a more extensive series 
of investigations, and the final discussion will be postponed to a later paper. 


(C. F. P.). 
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It was pointed out in (II.) that the new method is similar in principle 
to that employed in experiments with Wilson chambers operated in- 
magnetic fields. The average deviation in a track per unit length, @,. 
gives a measure of the quantity pB/Z for the corresponding particle, 
where p is its momentum, Ze its charge, and f the ratio of the velocity 
of the particle, v, to that of light, c. On the other hand, if the charge 
of the particle is known, the grain-density in its track gives a measure of 
its velocity—and thence the corresponding value of $. The two 
observations on the same track therefore allow both the mass and the 
velocity of the particle to be determined if its energy lies in a particular 
range of values ;_ viz., that its kinetic energy is less than 0-6 mc.°, where 
m is the rest-mass, and greater than ~15 MeV. The method has the 
advantage that it is possible to measure the specific ionization of a particle 
of the cosmic radiation, by observations on the “ grain density ” in its 
track in a photographic emulsion, with a much greater precision than by 
‘“ drop-counts ”’ from photographs taken with ‘ counter-controlled ~ 
Wilson chambers. Further, in favourable circumstances, it appears 
to be possible to determine the quantity pf/Z, by scattering measurements, 
with an accuracy comparable with that commonly achieved in experiments 
on the magnetic deflection of charged particles of the cosmic radiation. 

The importance of the new method lies in the fact that it greatly 
extends the field of application of the photographic method. Hitherto, 

using this method of detection, determinations of mass have been confined 
to particles of low energy which reach the end of their range in the 
emulsion. The new method removes this restriction and, within certain 
limits, allows us to determine the mass and energy of particles, irrespective 
of whether or not they escape from the emulsion. 

Whilst the early measurements, described in (II.), showed the 
potentialities of the method, the determination of the scattering parameter, 
z, was arduous, for it involved plotting the individual grains along a 
track from an image thrown on a screen by means of a projection micro- 
scope. This feature placed a limit on the field of application of the 
method ; for with the apparatus at our disposal, an extensive analysis, 
made desirable if observations of great statistical weight were to be 
obtained, would have occupied some years. It therefore became desirable 
to develop simplified methods; and we have found that the following, 
which employs microscopes of conventional design, has the advantage 
of greatly increasing the rate of measurement, without an appreciable 
loss in precision, in addition to extending the range of values of the energy 
of particles which can be determined. 


DETERMINATION OF THE SCATTERING PARAMETER. 


In order to determine the scattering parameter, «, for a given track, 
the plate is placed on the mechanical stage of the microscope, so that the 
track is approximately parallel to one of the lines of motion of the stage, 
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say the X-axis. The point, y,, at which the image of the track intersects. 
a hair-line in one of the eye-pieces, is determined by means of an eye-piece 
scale. The X-motion of the stage is then moved through a standard 
distance, say 100 1, and a reading taken of the new position, ys, of the 
Image of the track relative to the eye-piece graticule. The operations 
are repeated along the entire length of the track. 

The first differences, S;=y,—y,, S,=y,—y, ete. give a measure of 
the inclination 7, ij, ete. of the successive chords in the track—and the 
second differences, 8,—S,, S,—S,, etc. the deviations, D,, D,, etc.. 

n 


between successive chords. The mean value of |D|, =D, taken 


without regard to sign, gives a measure of the scattering parameter, 3, 
from which the quantity pf/Z can be evaluated. 

For the successful application of the method, it is necessary that any 
apparent deviations in the track due to mechanical imperfections of the 
stage motions shall be compared with the deviations due to scattering. 
it has been found that satisfactory performance in this respect is given 
by Type M4000 microscopes by Cooke, Troughton and Simms. These 
instruments can be fitted with a kinematic stage provided with micrometer 
adjustments on both the #- and y-motions. It appears to be difficult 
to make a simple test of the performance of a given stage without special 
optical devices, but observations onthe tracks of the particles with great 
momentum suggest that, when in good adjustment, any imperfections 
in the motion of the stage lead to errors which are not greater than those 
involved in reading the point of intersection of the image of the track 
and the eye-piece graticule. 

The results of a characteristic series of measurements on part of a 
track of length 8 mm. is shown in Table I. In this case, values of y have 
- been determined at successive points, 50 4 apart. The observations have 
been analysed taking a cell length, c, of 100, and the average value 
of ~ (100) the average angular deviation between successive chords 
of length 100, is found to be 0-105°. The individual values of D 
obtained by the above procedure, are not independent, but the use of 
overlapping cells tends to reduce the effects of reading errors. 

A similar analysis has been made using different values of c—200 p, 
300, etc. The mean values of « for different values of c are related by 
the equation «(c,)/%(c,)= V/(¢1)/+/(¢). The values of «(200 2), «(300 2), ete., 
can therefore be used to determine the value of «(100 1) to which they are 
equivalent, and the corresponding values are included in Table I. It 
will be seen that «(100 ) decreases with increasing cell size and reaches 
an approximately constant value. This is due to the fact that if the 
cell-length is too small, the apparent deviations in the track tend to be 
greater than the true values owing to reading errors, and the resulting 
value of « is too high. On the other hand, if the cell-size is unnecessarily 
large, the number of statistically independent observations is reduced, 
and the calculated values of « are subject to correspondingly greater: 
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statistical fluctuations. It is therefore desirable to employ the smallest 
possible cell-size which is consistent with the condition that the uncertainty 
in the value of % due to reading errors is small (Goldschmidt-Clermont 
et al. 1948). 


TABLE I. 


Characteristic results of the measurements of the deviations in the track of a 
particle of length 8 mm. Readings every 50 on first mm. of track. 
yrain-density in track 19-2 grains/50 np. c=100 p. 


& y ) D z y ) D 
mm dev Yn —Ynt 2 Sy Nar 2 mim. dev. Yn —Yn+2 Sn—Spt 2 
12-8 5:6 0-5 0-3 12-25 7-0 0-4 0-3 
12:75 5:7 0-6 0-6 12-2 TS: 0-2 0-2 
12:7 6-1 0-2 O-L 12-15 7:4 0-1 0-1 
12-65 6-3 0-0 —0-2 12-1 7-5 0-0 —0-3 
12-6 6-3 0-1 0-0 12-05 7-5 0-0 —():3 
12-55 6-3 0-2 0-1 12-0 75 0-3 0-3 
12°5 6-4 Dell —0-2 11-95 7-5 0-3 0-3 
12-45 6-5 O-1 —0-3 11:9 7:8 0-0 —(-1 
12-4 6-5 0-3 —0-2 11-85 7:8 0-0 —()-] 
12-35 6-6 0-4 0-0 11-8 7:8 O-l 0-0 
12°3 6-8 0-5 0-3 


Values of a(100 4) for different values of c. 


| 
a(100 p)° 


5 fame 
microns 


50 0-141 | 0-191 
100 0-22 0-105 
150 0-285 | 0-074 
200 0-35 0-060 
300 0-54 0-050 
400 0-78 0-047 
600 1-56 0-051 

1000 3°74 0-057 
1400 6-28 0-058 
2000 11-0 0-058 


For small values of ¢, it is found that, owing to reading errors, values 
of D tend to become constant, and independent of .¢ It will be convenient 
to follow previous authors and refer to this as the “ noise-level ” (Dy 
1? € re . a a . 
lhe actual values of D,, depend on a number of factors such as the grain- 
density in the track, the quality and adjustment of the optical system 


of the microscope, and the individual skill of the observer, but they are 
-of the order of 0-2 p. 
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~ “As the cell-size, c, is increased from these low levels, the value of (D) 


rises above the noise level (D,,). The corresponding values of «(100 y). 
then rapidly approach a constant value apart from statistical 
fluctuations. These effects may be summarized 1 by saying that the errors. 
of measurement become unimportant if D> 4D,. The minimum value 
of ¢ which gives this condition is adopted in the deter mination of «(100 4): 

It i is clear that if the value of « is to have a given statistical weight, the 
minimum length of track accepted for measurement will pend: on the 
value of «(100 «), long tracks being necessary for particles of great energy. 
Thus if Lis the length of the track, the number of statistically independent 
values of D is (L/c)—1. It can be shown that the root-mean-square of the 
deviations of the individual values of D, from the true mean value corres- 
ponding to the momentum of the particle, is equal to ~0-5D. Approximate 
values of the minimum length of track required to ensure that the probable 
error in the determination of D is less than 25 per cent, are included, 
for different values of «(100 2), in Table IT. 


TABLE II. 


Values of the minimum length of a track which allows the scattering parameter, 
_ a(100n), to be determined with a probable error of less than 25 per cent, 
_ in the absence of distortion. 


Length of track a(100p) 
mm. degrees 

3 0-017 

5 0-008 

10 0-003 

20 ; 0-001 


Errects oF DISTORTION: OF THE EMULSION. 


When an exposed plate is processed, there is a large increase in the 
volume of the emulsion. If there is good adhesion between the gelatine 
and the glass, the swelling in the central regions of the plate corresponds 
approximately to a uniform linear transformation in the Z co-ordinates 
for any point ; in the direction, that is, perpendicular to the glass. After 
fixation and washing, the drying of the plate produces the well-known 
shrinkage associated with loss of water and theremoval of the silver salt from 
the emulsion during fixation. Since we are concerned only with deviations 
in the plane projection of a track, such a shrinkage will not, in principle, 
affect the measurements. In practice, however—and especially with 
thick emulsions, and in regions near the edges of the plate—some distortion 
of the emulsion occurs, so that a strictly rectilinear track will, in general, 
appear to be curved. 

‘The distortion of the emulsion appears to arise, at least, in part, as a 
result of stresses released during the absorption of water, and because the 
top surface and sides are unconstrained. An emulsion, of which the cross- 
section in a plane perpendicular to the glass is rectangular when the 
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emulsion is dry, becomes distorted, therefore, when placed in the processing 
solutions. Further, the distortion may persist or be increased when the 
plate is dried. The difficulties can be greatly reduced by using large 
plates and confining the observations to the central areas ; but even with 
these precautions, deviations due to distortion can give rise to serious 
errors in the determination of the energy of the fastest particles. 
(a(100 4) <0-01° per 100). In the worst cases, corrections may be 
necessary with tracks for which the true value of «(100 ) is as great 
as 0-05°. 

Evidence for the existence of distortion can be obtained in two ways : 
first, it is sometimes found that the deviations, D, instead of being 
approximately equally distributed between positive and negative values, 
are all of one sign. In such cases, an inspection of more steeply dipping 
tracks, of which the projection is approximately parallel to that of the 
track under examination, can show whether distortion is indeed present. 

Observations show that the distortion is commonly of such a character 
that the deviation which it produces in a strictly rectilinear track can be 
written as 4=Féz, (1), where 4 is the change indirection in the track 
in passing from a depth z to z+6z in the emulsion, and F is a constant 
for a particular direction of the projection of the track in a given region 
of the plate. From observations on a more steeply dipping track, the 
constant F can be determined and appropriate corrections then made to 
the track under examination. If these corrections change the values 
of the deviations attributed to scattering so that they are now (apart 
from differences associated with statistical fluctuations) symmetrically 
distributed about the value zero, the corrected values are accepted as 
adequate for the determination of the true value of D. 

A better method, and the one most commonly adopted, is the following : 
distortion of the type represented by equation (1) has the result that the 
apparent deviations in the track are due to the joint effect of the deviations 
due to scattering, D, and those due to distortion, ¢. The apparent 
deviations, D,, may therefore be written as D,=D-+¢, where ¢ is a 
constant term*. Ifthe distortion is severe, so that ¢ is large, the values 
of D, will be all of one sign in spite of the fact that values of D are equally 
positive and negative apart from statistical fluctuations. 

An approximate value of ¢ may be obtained from the quantity 
2D,/n~, where 2D, represents the algebraic sum of the individual 
values of D,. With this value of 4, the individual values of D can be 
found and hence an approximate value of |D|. The value of |D] thus 
obtained is insensitive to the actual value of 4 employed provided the 
error in ¢ is less than | D|/2. 

Table III. shows the results of measurements on a track which shows 
evidence of severe distortion. In addition to the correction ¢ calculated 
by the method described above, a number of arbitrarily chosen values 


*It is a sufficient approximation to assume that the track “dips” at a 
constant inclination to the plane of the emulsion. 
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have been employed, and the results summarized in Table III. illustrate 
the fact that the final result is insensitive to the precise magnitude of the 
correction term. 

TasBie ITI. 


Observation on track showing evidence of distortion. 


: . D* 
iG ) De Geek 
mm divs Yn—Ynt+e2 Si =S7h 2 (Ome 

9-1 10-0 —0-5 5-0 1-0 
8-9 9:5 1-0 7-0 3-0 
8-7 9-5 4-5 6-0 2-0 
8:5 10-5 8-0 4-5 0-5 
8:3 14-0 10-25 5-0 1-0 
8-1 18-5 12:5 5°25 1:25 
7:9 24-25 15-25 4-0 0-0 
7:1 31-0 17°75 2-0 —2-0 
(5) 39-5 19-25 1-25 —2-75 
lees 48°75 19-25 2-15 —]-25 
ae] 58-75 20-5 3°25 —0Q-75 
6-9 68-5 22-5 2-0 —2-0 
6:7 79-25 23°75 
6-5 91-0 24-5 
6:3 103-0 
6-1 115°5 


d=ZDin=4; U2D*/n=1-46 


Values of «*(100 «) for different values of ¢. 
$ | 3D—4) | ™(l00p) 


0 48-0 0-12 
I 36-0 0-09 
2 25-5 0-064 
2) 20-0 0-050 
4 17-5 0-044 
i) 18-5 0-046 
6 26-0 0-065 
iT 36-0 0-090 
8 48-0 0-12 


LIMITATIONS OF THE METHOD. 

The minimum value of D which can be regarded as significant, and 
+herefore the maximum value of the energy of.a particle which can be 
measured, depends on the length of the track available for measurement, 
the magnitude of the “ noise-level ”, and the degree of distortion of the 
gelatine. The maximum measurable values of the energy corresponding 
to tracks of different length given in Table II., correspond to conditions 
in which distortion is absent. The minimum value of o(100~) actually 


G2 
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observed in the present experiment is 0-001°, and we have no evidence 
+o show whether the theoretical figure of «(100 4)=0-0005° for a track of 
length 2cm., can in fact be attained in our experimental conditions. 
Such an observation, to be significant, would demand that the motion 
of the stage departs from being strictly rectilinear by an amount less 
than 0-4 1 in traversing a distance of 2cm. That a performance of this 
order of precision is obtained for shorter lengths of travel is proved by the 
observation of long tracks with a value of «(100 4) =0-001°. 

If distortion exists, then the minimum value of D which can be accepted 
as significant for a track of a given length is raised above the values. 
given in Table II. By confining the observations to tracks in the central 
regions of the plates, however, and with the convention that the distortion 
correction ¢ is less than 2D, the number of tracks which are rejected when, 
but for distortion they would have given significant results, is very small. 


CALIBRATION. 

It has been shown (E.J. Williams 1939, 1940; Rossiand Greisen 1941) that 
if the track of a particle is divided into cells of length c, the average 
angular deviation between successive chords «, is related to the velocity »,. 
of the particle by the equation 


i Huo{(,fiogtas)h, a) 


where Ze is the charge on ne particle, and m, its rest-mass. The constant 
in equation (2) depends on the composition of the emulsion, and on the 
chosen length of cell, c. It is a sufficient approximation to neglect 
changes in the value of the logarithmic term in (2). 
k,Z Z 

—; (1— P= he — oe a SE) 


Mor” pe 


o= 


where mp is the rest-mass of the particle and p its momentum. 

In order to determine the relation between the grain-density, g, in a 
track, and the rate of loss of energy of the particle producing it, observations. 
have been made on the tracks of mesons, protons and deuterons of great 

range which stop in the emulsion. The results are shown in fig. 1. This 

curve, and the range-energy relation for the same particles, would allow 
us to determine the grain-density in the track of any particle of charge e, 
as a function of its velocity. | 

Direct determinations of the range-energy relation for protons in 
“ electron-sensitive ”’ emulsions have not yet been published. It is 
known, however, that in the Ilford C2 emulsions, the range R, in 
microns, of a proton with an energy E MeV., is given by the relation 

R=CE’ 


and that dH/dR=kz*p" (4). Assume, therefore, that in the G5 emulsion: 
arelationship of a similar ae exists. The corresponding values of k and. 
n can be determined from the following observations : 
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Relation between grain-density and rate of loss of energy along the trajectory. 
Ilford G.5 emulsion. 


in kev/p. 


(1) The average range in this emulsion of p-mesons, emitted with a 
kinetic energy of 4:1 MeV. from a 7-meson at rest, is 595+ 10 pe. 
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(2) From fig. 1, we determine the residual range of a proton which has. 
an initial ionization, as determined from the grain density in its track, 
equal to that of energetic “ primary ” «-particles with a specific ionization 
equal to four times the minimum value for charge je]. 

From these two observations it is found that 


dE 058722 


aR ras Kev. per micron. eT) 


This value of » is identical with that for the C2 emulsion. 

Using this relation, together with the variation of the grain-density 
along the track of protons shown in fig. 1, we can determine the grain- 
density in the track as a function of the rate of loss of energy of the particle; 
see fig. 2. It will be seen that the relationship is almost linear for values. 
of the rate of loss of energy up to four times the minimum value. From 
equation (5) and fig. 2, the relation between the grain-density in the track 
of an energetic proton, and the corresponding value of the scattering 
parameter «%, can be determined. 


.XPERIMENTAL RESULTS. 


The methods employed in the previous paragraphs have been applied 
in a study of 350 tracks, associated with “stars ’’, and chosen only 
because their length in the emulsion is greater than 3000. The 
observations were made in plates exposed at high altitudes as described 
in Part II. of the present series of papers. The results are represented 
in fig. 3 where the value of the scattering parameter, x, for a given track, 
is plotted against the corresponding value of the grain-density, g, both 
quantities being on a logarithmic scale. 

The lines in fig. 3, marked z, P, D, and T are of the form to be expected 
for particles of charge e as determined by the methods of the preceding 
paragraphs. They are, however, displaced by about 7 per cent from the 
calculated curves for the corresponding particles along the axis of the 
scattering parameter. Their actual positions have been determined 
by the following method : 

The calculated curve for protons has been displaced until it appears 
to give the best fit with the points corresponding to the most prominent 
group of particles. These: are certainly due to protons. The other 
curves, with the same form as the first, are then drawn to correspond to 
particles of mass 286 m,, for 7-mesons; and with appropriate values. 
of the mass for deuterons and tritons. The corresponding curves for 
z-particles, lithium nuclei, ete. can also be deduced from the curve for 
protons using the information given in fig. 2. The difference of only 
7 per cent between the calculated curve for protons and the observed 
distribution for these particles is not unexpected in view of the uncertainty 
in some of the quantities on which the calculations are based, but it. 
remains to be explained. 
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In plotting the results shown in fig. 3, primary and secondary particles 


have been distinguished by the criteria described in (I.); (Brown et al. 


Log, Grain-density, grams per 50 p. 


1949). In the case of events in which many shower particles are emitted, 
a particle which produces a track with a direction near the axis of the cone 


Fig. 3. 
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Relation between the grain-density and a(100), for the tracks of 350 particles. 
“Primary ” particles are shown thus @ ; and “‘ secondary ”’ particles thus ©. 
The line P has the same form as the calculated curve for protons, but is 
displaced to the left, in the direction of the a axis, by 7 percent. The lines 
“a? T)” and“ T’”’ are drawn in the correct relation to line P to correspond 
to particles of mass 286 m,, 3674m, and55llm,. “A”, “B” and“ C” are 
lines of constant momentum, and correspond. to 500, 1000 and 2000 MeV./c, 
respectively. Energy scales, in MeV., for z-particles and protons are 
included at the top of the figure. It will be seen that the a-particles of 
great energy are all “ primaries’’. The two points corresponding to particles 
of greater mass are due to a “ primary ’’ nitrogen and the “ secondary ” 
lithium nucleus occurring in the same event. 
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embracing the secondary particles, but on the opposite side of the star 
from the “shower” particles, is regarded as the primary. In the 
case of events of lower multiplicity, and especially in observations :att 
great altitudes where the directions of motion of the primary particles 
are more widely dispersed than at 11,000 feet, the identification is subject 
to some uncertainty, and must be-treated with reserve. The greater 
values of the energy of the primary, as compared with the secondary 
particles (see fig. 3) strongly suggests, however, that the errors in 
identification are not serious. 

It will be seen, from an inspection of fiz. 3, that the primary «-particles 
produce disintegrations with a frequency of the same order or magnitude 
as the protons in the same range of energies. On the other hand, compared 
with protons and deuterons, «-particles of great energy are emitted from 
nuclear disintegrations very infrequently. 


Fig. 4. 
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Distribution in mass of the particles of charge [e|, corresponding 
to the measurements represented in fig. 3. 


_ The displacement from the proton curve, in a direction parallel to: the 
a-axis, of any observed point in fig. 3 gives a measure of the mass of the 
corresponding particle. The distribution in the values thus obtained— 
based on measurements of all the tracks represented in fig. 3 which have 
values of the grain-density lying in the interval from 15 to 98 grains 
per 50 p— is shown in fig. 4. It will be seen that in addition to the 
prominent peak due to protons, there are groups of smaller intensity, 
unresolved from the protons, which may be attributed to deuterons and 
_tritons ; and in addition, a well-defined group due to mesons of which the 


mean mass is 2834-14 m,. 
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_ it may be emphasized that if, as an alternative procedure, the dis- 
placement of the individual points from the theoretical curve for protons 
had been measured, the resulting distribution in the values of the mass 
of the particles would have been identical in form with that shown in 
fig. 4, but all the values would have been lowered by 7 per cent. It is 
reasonable, however, to deduce the mass of the mesons by a direct 
comparison with that of protons, and the result obtained is identical with 


that given above.. The majority of the mesons may therefore be identified 
-as 7-particles 


Number of tracks. 


Grain-density (per 100 y). 
Pistribution in the values of the grain-density of the tracks of particles of 
charge Je]. The lower distribution: is for all particles ; and the upper, for 
those which are identified as z-mesons. 


The observed distribution in the values of the grain-density of all the 
tracks of particles of charge ¢ represented in fig. 3 is shown in fig. 5. This 
curve may be compared with the results of similar measurements made on 
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plates exposed on Jungfraujochand givenin (I.), fig. 1. The two distributions. 
are similar in form, but there are relatively more particles of low grain- 
density in the present observations. This corresponds to the greater 
average energy of the nuclear explosions at the greater altitude, see (Ee 

A study of fig. 3 shows that most of the particles producing tracks with 
a grain-density less than 14 grains per 50 p are either mesons or ~ primary ”’ 
particles. The nature of the secondary particles for which «(100 ;.) <0-02° 
cannot be determined by our present methods. It is possible, however, 
to extrapolate the observed distributions in energy of protons and 7-mesons, 
in the range in which they can be distinguished, to the high energy region,. 
and thus to estimate the proportion of the two types. Making suitable 
corrections for the “‘ escape” of particles from the emulsion, the relative 
frequency of occurrence of the different types can be estimated, and the 
results are shown in Table IV. 


TABLE LV. 


Number of particles of different types, per “ star’’, with a range greater than 
3 mm., observed at an altitude of 70,000 feet. 
m-particles 0-6, Protons 1-4, 
Deuterons 05, Tritons 0:2, 
a-particles 0-02. ‘ 


Average number of tracks per “ star’ =7-4. 


CONCLUSIONS. 


The following conclusions can be drawn from the results of the present 
experiments : 


(1) Of the mesons emerging with a kinetic energy less than 150 MeY. 
from nuclear disintegrations produced by protons and «-particles of great 
energy, the great majority are z-particles. If y-particles are sometimes. 
directly created, they are present in numbers less than 2 per cent of the 
7-particles. 

(2) If electrons, or other charged particles of small rest-mass, are 
sometimes emitted from nuclear explosions, with an energy less than 
150 MeV., they occur with a frequency less than 2 per cent of the z-particles.. 
Most of the electrons of the soft component of the cosmic radiation have 
an energy less than 200 MeV. The observations therefore prove that 
most of them are of secondary origin. 

(3) The observations give no support for the view that mesons of 
intermediate mass, 7-mesons wi.h a lifetime greater than 10°™ sec., are 
ejected with an energy less than ~300 MeV., during nuclear explosions ; 
and they prove that if such particles, with a mass between 350 and 
1200 m,, are sometimes created in these processes, their number is less 
than 2 per cent of that of the z-mesons and less than | per cent of the 
protons. This conclusion is complementary to that drawn recently by 
Franzinetti (1950) from determinations by magnetic deflection experiments 
of the mass of slow charged particles produced by cosmic radiation. 
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(4) No final conclusion can be drawn regarding the nature of the particles, 
less massive than the proton, which produce tracks with the minimum 
value of the grain-density. In view of the observations on particles of 
lower energy, it is reasonable to suppose that the faster particles are of 
the same type. In what follows they will therefore be regarded as. 
7-particles. 

(5) At least 90 per cent of the * shower” particles emerging from nuclear 
explosions are 7-mesons, a “‘ shower” particle being defined as one with 
a grain-density less than 1-5 g,,;,, (see I.). Consider, however, all particles 
im a given range of momentum—say between 200 and 700 MeV./c—a 
practice commonly adopted in experiments with Wilson chambers 
operated in magnetic fields. Among the particles in this momentum 
range, the 7-mesons are accompanied by large numbers of protons. If 
the a-mesons are equally distributed between positive and negative: 
_ particles, the “ positive excess” should be of the order of six, a value 
consistent. with that observed in experiments with expansion chambers. 
(Rochester 1949). 

(6) An extension of the observations will allow a detailed analysis of 
the nuclear processes produced by protons and «-particles in the range: 
of energies above 10° eV., including the distribution in energy of the 
7-mesons near their points of production, and the relation between the energy 
ofa 7-particle and its direction of emission with respect to the * primary ” 
particle. Such investigations are being undertaken in a collaboration with 
U. Camerini, W. O. Lock and H. Muirhead; and by J. H. Davies and 
C. Franzinetti. 

(7) The observed homogeneity in the values of the mass of the ejected 
mesons with energy less than 150 MeV. allows us to exclude certain 
assumptions regarding the origin of the particles. It is known, from the 
experiments on the artificial creation of 7-mesons made at Berkeley 
(Gardner and Lattes 1949), that the threshold for the process occurs at 
such an energy of the bombarding protons or z-particles that the creation 
of a particle of mass ~280 m, is just possible from the point of view 
of the conservation laws. It follows that the 7-mesons must be 
directly created ; that they cannot be regarded as the products of decay 
of mesons of appreciably greater mass, and very shoré lifetime. 

In the nuclear transmutations with which we are concerned, limitations. 
of energy do not apply. In principle, therefore, the possibility cannot 
be excluded that the 7-mesons occurring in showers are produced by the 
decay of short-lived mesons of greater mass which are the primary 
products of the nucleon-nucleon collisions. The observed homogeneity 
in the values of the mass of the ejected mesons proves, however, that any 
postulated unstable parent of the 7-particles, if it exists, must be very 
short-lived; for if it survived, on the average, over a length of path of 
say 2 mm., the observed values of the mass would represent averages 
between the masses of the parent and daughter particles. As a result, 
we should expect to observe a dispersed set of values instead of the compact 
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group actually found. The present observations therefore prove that 
any more massive parent mesons have lifetimes less than 107% see. 
It is, of course, reasonable to suppose that the z-mesons are directly 
produced in the nucleon-nucleon collisions. 
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ABSTRACT. 

The ati is considered of obtaining a divergence free S-matrix 
for the interactions of mesons with nucleons along the lines of mass and 
charge renormalization, which Dyson has so successfully applied to 
quantum electrodynamics. It has been shown that this can be done 
for the vector interaction of neutral vector mesons. The procedure is 
here extended to the scalar interactions of scalar and pseudo-scalar 
mesons but proves to be quite inadequate to remove all the divergences. 
im the other cases. 


§ 1. IyrRopvucTIoNn. 


By applying Feynman's (1949 a, b) graphical methods for the calculation 
of matrix elements to relativistic quantum electrodynamics in the 
Schwinger (1948) form, Dyson (1949a, b)+ has shown how a finite 
S-matrix can be obtained for a wide class of phenomena. In Feynman’s. 
theory, any real process taking place with specific virtual processes can be 
represented by a graph from which the integral for the corresponding 
matrix element can be written down directly, each line and each vertex 
corresponding to a single factor in the integrand. An essential feature 


‘of Dyson’s development is the proof that only a limited number of 


different types of divergence can occur in quantum electrodynamics. 
These are called primitive divergents. When a primitive divergent graph 
occurs as part of a larger graph its effect is to complicate a single line or 
a single vertex of a simpler graph, and consequently to modify a single 
factor of the integrand corresponding to the simpler graph. Ifthe infinite 
parts of all such modifications of a graph are dropped in a well-defined 
way, a finite S-matrix is obtained for all phenomena which do not involve 
bound states. The mass and charge of a free particle in the vacuum are 
both given by the theory as infinite renormalizations of the mass and 
charge of an ideal ‘ bare’ particle. Dyson has shown that his prescription 


for dropping the infinite terms is equivalent to equating the infinite 


theoretical values of mass and charge to the finite experimental values. 
The extension of this method to the various meson-nucleon interactions. 

is here considered. The number and degree of the possible primitive 

divergents in any theory are closely related to the degree of the derivatives 


* Communicated by Dr. N. Kemmer. | 
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.of 4-functions which occur in the commutation relations of the Hamiltonian 
under consideration. The meson interactions which in this respect most 
closely resemble electrodynamics are the scalar interactions of scalar and 
pseudo-scalar mesons with nucleons in which no derivatives occur from 
meson variables. The possible types of primitive divergents are here 
identical with those of electrodynamics. In electrodynamics one of these 
—the scattering of light by light—is finite by gauge invariance, but the 
scattering of mesons by mesons gives rise to genuine divergences. However, 
it will be shown that these can be regarded as a special type of self-energy 
between pairs of mesons and can be cancelled by adding the appropriate 
.6-type interaction between mesons to the Lagrangian. A finite S-matrix 
will thus be obtained. In the remaining meson-nucleon interactions 
double derivatives arise from the commutation relations of the meson 
variables with the result that there is no limit to the number and degree 
‘of the possible types of primitive divergents. Thus the renormalization 
-of mass and charge is shown to be quite inadequate for the removal of 
all divergences in these cases. The only exception to this is the vector 
interaction of neutral vector mesons with nucleons, which has already 
been dealt with elsewhere (Matthews 1949 b). . 

The proofs of this paper do not confirm exac'ly the conclusions of Case 
(1949) which were generalizations from the resul_s of particular calcula‘ions. 
Some of the infinities considered here only occur in comparatively high 
order correciions and were thus not detecied. 


§2. THE S-MATRIX. 
The interaction of two fields can be expressed by a- generalized 
Schrédinger equation in the interaction representation, 


ihes P[o]/do(z)=H(d(a)) Plo], 2 8 ee 


where W[o] is the wave function on a general space-like surface, o, and 
f(x) is a typical field variable satisfying the equation of motion and the 
commutation relations of its free field with mechanical mass. The 
Hamiltonian in (1) has matrix elements for transitions which cannot 
occur with conservation of energy and momentum, and it is convenient 
to remove these by a Schwinger (1948) transformation 


Y+>exp{—iS(d)}¥. JT Be fo, ee SO 
‘The new equation is 
ihe(3/dc)(exp{—iS(¢)}¥)=H(4) exp{—iS(p)}¥ 2 2. (3) 


‘To second order in the coupling constant the new Hamiltonian has matrix 
elements for real processes and also some infinite elements for the self- 
energies of the interacting particles. If the infinite self-energy parts 
of the Hamiltonian in (3) are H,(¢), it has been shown for a fermion 
field y% and a boson field ¢,, (Matthews 1949 ¢), that 


(4, b)=hediegipab+ hSx2 eat Se aie ee eres 
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where xk, and « are the Escrerovall Compton wavelengths for the fermions 
and the bosons respectively. To obtain an equation from which these 
terms have been removed make the transformation 


| ee ee coe Sys Re NS PRS) 
where 
hcedU/do=H_,(d)U. ery Pee KG) 
It then follows that 
UG Wah Bet om & BAG) aod Oa eee UB 
where 
ae = ee ee aa, : (8) . 


The equation which Dyson makes the starting-point of his ee 
obtained by reversing the order of these transformations. The 
transformation (5) on equation (1) leads to 


- ihe(S¥ ./5c)-+-ihcU-1(8U/8c) ¥ ,=H(¢,)¥,. 
y (7), 
ihed¥ 0 {H($,)—H(.)}¥, | 
Se) Oe Mo iare ete es (0) 


which is the Schrédinger equation in the interaction representation in a 
. form which is free from self-energy effects up to the second order in the 
coupling constant. Since ¢ satisfies the free field equation with mechanical 
mass, it follows from (4), (6) and (8) that ¢, satisfies the free field equation 
with renormalized mass. This is given infinite by the theory, but since 
it is the mass of the observed free particle, it is equated to the finite 
experimental mass. An equation like (3) with the self-energy explicitly 
subtracted could now be obtained from (9) by performing the 
transformation (2) with ¢ replaced by ¢,. 

Dyson’s theory starts with an equation of the form (9), but the constants 
Ok, and 8x? are chosen as power series in the coupling constant so that 
they will cancel all the self-energy effects from the rest of the Hamiltonian 
up to any order. If Y% is the state vector of a system on a surface in 
the remote past and YW is the same system, after its particles have 
interacted, in the distant future, then 


eS cn | Ure ae pe ee ted (0) 


where in Dyson’s notation (I. §5; II. § 2) 


S(0)= = £(E) x : fo ae. Fas dae POEL (2), Ha), oo Hf) (11) 


‘S( co) is the Heisenberg S-matrix, and, by developing Feynman's (1949 a, b) 
methods, Dyson has given comparatively simple rules for the calculation 
of its elements (I. §§ 6, 7; II. § 2). Essentially, if H, contains factors 


x, ib and ¢, where 
Poteet isu A( arm ramet deta <CL2) 
and 


{b(x), (y)}= —ty,,0/0x,,—Ky)A(a@—y). 6. . (13) 
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the matrix element for a process in which a virtual particle passes between 
a and y, contains in the integrand a factor 


(P(b(«), 6y))o=HicAp(e—y) . . . . « (4) 


if the particle is a boson, and 


(Pah (x), H(Y))) o= — 40%, Y)(%4.8/ 0%, +) Ap (@—Y) = — 39(@, Y)Sp(%@—Y) 


if it is a fermion. If, in place of 4d, H, contains a factor 0¢/dx,, then it 
will also contain a term depending explicitly on the normal to the surface 

co. It has already been shown (Matthews 1949 a) that to extend Dyson’s 
‘ rules to this case the normal-dependent terms in H, can be neglected 
and then effectively 


(P(24/2x,, 94/8y,)) =e 22/8x,8y,)Ap(a—y). . - ~ (15) 
A similar result holds if the derivatives occur in the commutation 
relations themselves. Thus for vector mesons 


[$,(@), 6,(y)l=the(8,,+ (1/K?)0?/dx,0y,)A(a—y), . . (16) 


and to calculate terms of (11) for a Hamiltonian containing a factor ¢,, 
the normal-dependent terms can be neglected and effectively 


(P(,(%); b4(Y))) o=BAC(Spy + (1/2) 02/Ox,0Y,)Ap(w—y). (17) 


The infinities in 8( 00) associated with mass renormalizations have already, 
in principle, been removed. There remain other infinities which in 
electrodynamics can be removed entirely by renormalization of the charge. 
Before considering the possible application of this method to meson 
interactions, it will be necessary to discuss in more detail Dyson’s analysis 
of the possible types of primitive divergents, because it is at this point 
that the difficulties of extending the theory arise. 


§ 3. PRIMITIVE DIVERGENTS. 


Consider an interaction between fermions and bosons in which no. 
derivatives occur either in the interaction or in the commutation relations 
of the bosons. Every element of S( 00) is associated with certain graphs 
G, the process under consideration determining the number of external 
fermion lines, K;, and the number of external boson lines, E,. Consider 
any such graph, G, with » vertices and F internal lines of which Fy are 
fermion lines and I, are boson lines. Now 


yi mh. fexp{—ip,%,) 
y(2)= Gas [EE Oy, dsp 
The integrals over the space variables in the elements of S( 00) are trivial, 
giving conservation of energy and momentum at each vertex. Thus 
the energy-momentum associated with each line in G can be expressed 
in terms of (F—n--1) basic vectors p', and the corresponding contribution 
to S(0o) is an integral, M, over these variables. For considerations of 


(18) 
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convergence, each boson line contributes to the integrand of M a factor 
1/((Xa;p")?+-«?) and each fermion line a factor Xa;pi./((Xa;p')?-+-«2), where 
a;p' is the energy-momentum associated with the line. The numerator 
of the fermion factor comes from the derivative of the A-function in the 
fermion commutation relation (13). Many of the coefficients a, will be 
zero and, if for a certain line only a, is non-zero, the line will be referred 
to.asa basic line. If an internal line ina graph is replaced by two external 
lines, the line is said to be opened. 

The differentials d‘p,, can be expressed in polar variables. The integrals 
over the angle variables are finite and the integral over the radius is p3dp. 
Define the convergence function 


K=2F—F,—4(F—n+1). 2 2. 2 2. (19) 


If there is only one basic p‘, —(K-+1) is the power of p in the integrand 
of M expressed in polar variables. The condition for convergence is 
in this case 


elmer AP Tero obo. (20) 


If there are several basic variables, define K, to be K evaluated for the 
graph obtained by opening all the basic lines except p’, that is for the 
integral M with all variables except p‘ held constant. Similarly define K ;, 
to be K evaluated for the graph with all basic lines except pé and p? 
opened. A sufficient condition for M to be conditionally convergent is 
that 

Se Wen Kee say ae Bot a) grater gr(21) 
for all 7,7,...=1,2,...,(F—n+1). The graph is absolutely con- 
vergent if (21) holds for all possible choices of the basic variables p’. 
This condition is not necessary for convergence because it may happen 
that the coefficients of the higher powers of p in M vanish identically. 

A graph is said to be primitive if the graphs obtained from it by opening 
any line are convergent. Thus a primitive graph certainly satisfies all 
the conditions (21) except possibly the last one. If this is also satisfied, 
the graph is itself convergent. The only possible primitive divergent 
graphs are those for which the last condition is not satisfied. That is if 


SS Ee ey eT mea 


If G has n, vertices with just two fermion lines incident, and », vertices 
with just two boson lines incident, then 


IK —3n Me Np i; E;,, ‘ C co . C 5 (23) 

and ! 
Fp=n—% —$H,. eet cen) (24 p 

Thus , 
K=(3/2)H, t- H,4 Ny + 2m —4, Sle” Me aie Sanaa 6159 


and the only possible primitive divergents are (i) H;-=2, E,=0; (ii) H,=2, 
E,=1: (iii) E,=0, E,=2; (iv) E,=0, E,=4, and two other graphs which 
contain odd fermion loops which cancel with the similar graphs obtained. 
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by reversing the directions of the lines in the loop. The graphs (i), (ii) and 
(iii) are self-energy and vertex parts which can be regarded as modifications 
of single lines or single vertices. The graphs (iv) represent the scattering 
of bosons by bosons and will be called ‘‘ square * parts. The simplest, 
containing four free boson and a loop of four fermion lines, will be called 
a simple square part, Sj. 

Now any divergent graph must be a divergent primitive or contain 
divergent primitive parts and a convergent graph can only be made 
divergent by the insertion of primitive divergent parts. Dyson showed, 
(IL. § 6), by general arguments of gauge invariance that for electrodynamics 
no new infinities arise from the insertion of square parts. These arguments 
are of course not applicable to mesons, and by Dyson’s methods it can be 
shown by direct calculation that the scattering of scalar or pseudo-scalar 
mesons by themselves is not finite like the scattering of light by light. 
Thus for meson interactions a special consideration of square parts is 
necessary. 

For a single free meson in the vacuum the self-energy arises from the 
virtual creation of nucleon pairs. When there are two free mesons these 
can scatter each other by virtual nucleon pair creation, (square graph). 
The matrix element for this process is infinite, but a finite part can be 
separated off if, in the spirit of Dyson’s work, we subtract from it that 
part which is obtained by setting the energy-momenta of the free mesons 
equal to zero. This subtracted part can be regarded as a special type of 
self-energy, which will have the form hcdA¢*, where 6A is an infinite 
constant. In the next section, it will be shown how a finite S-matrix 
can be obtained up to any order for certain meson interactions by the 
addition of a term of this type to the Lagrangian for the interaction. 

Note that the 5A effect, unlike the 8x”, is essentially an interaction 
between mesons which will not take place for a single free meson. The 
term is thus properly included in the interaction Lagrangian, when it 
will occur as an additional term in the Hamiltonian of (1) and, after the 
transformation (5), in the Hamiltonian H,(¢,) of (9). It would not be 
a reasonable procedure to start with the usual Lagrangian and to add 
the term to the self-energy Hamiltonian, H,, of (4). It would then 
appear both in H,(¢,) of (9) and in the free field Hamiltonian determining 
the motion of ¢,. But the equation of motion of ¢, should be the wave 
equation for a single free meson which, as pointed out above, should not 
depend on 6A. 


§ 4. ScaLAR Mgson SCALAR INTERACTIONS. 


The only interactions, in addition to electrodynamics, which satisfy 
the condition at the beginning of the last section are the scalar interactions 
of scalar and pseudo-scalar mesons with nucleons. If the term fhcdAd4 
is added to the Lagrangian in the pseudo-scalar case, the Hamiltonian 
H,(a) of (9) is 


H (ej iff (ae)yyguf(ar oh x) — icSaqys(x)yb( 22) — bdx*d?(x)—hedAd*(x). (26) 
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‘This Hamiltonian will give rise to graphs in which the most important 
vertices will have two nucleon and one meson line incident, at which 
the first term of (26) is operating. There will also be vertices at which 
either two nucleon, two meson or four meson lines are incident. At these 
vertices the terms dx, d« and 5A respectively will be operating. The rules 
for obtaining the matrix elements from the graphs are the same as those 
given by Dyson, (IT. §§ 2, 3), with the following substitutions. A,(2) is 
replaced by 4(x), ¢ is replaced by —f, y, is replaced by y; and D,(a) is 
replaced by 4,(wv). Also if n, is the number of vertices at which, 5A is 
operating, in place of IT. Eq. (12), terms must be multiplied by 
em ar eres POy 8 a Se VPS tot (2'7) 
‘This term makes allowance for the changes of sign which are due to the 
rearrangement of the and # factors to bring associated pairs adjacent, 
(1. §§ 6, 7). As in electrodynamics, functions X(W), I7(W’) and A(V) 
can be defined, (II. Eqs. (27), (39), (32)) through which can be determined 
the alteration in the integrand due to the insertion of a nucleon self-energy 
part, W, a meson self-energy part, W’, or a vertex part, V. If W consists 
of a single vertex, 
We theme OOK ne es) eh ie ae. (28) 
If W’ is a single vertex, 
NNO MG eamlOKe ob ee ek tt (29) 
The different signs of the terms in (28) and (29) are due to the factor (27). 
It is also necessary to consider explicitly the integrals from square graphs, 
S, which will be of the form 
Ap (t)4_()4,(2)4,(4-+-?—#)K(S, #, ?, ¢). 
In general, K(S) is itself a function of 4,, S, and y;. In the special case 
when § is a single vertex, 
See tr) =etON sy cA eget) cSt 1380) 
In Dyson’s terminology, a “ proper ”’ vertex or self-energy part is one 
which cannot be split into two parts joined by a single line. An 
‘irreducible ’’ graph is one from which all self-energy and vertex parts 
have been removed. It is also convenient to consider irreducible self- 
energy and vertex parts, which are such parts that are proper and from 
which all internal self-energy and vertex parts have been removed. 
We now extend these definitions. A ‘‘ quasi-proper”’ square part is a 
‘square part which cannot be split into a square part joined to a self-energy 
part by asingle line. <A ‘‘ proper ” square part is a quasi-proper square 
part which cannot be split into two square parts joined by two meson 
lines. ‘‘ Totally irreducible ’’ graphs are irreducible graphs which contain 
only simple square parts. X(t) is the sum of 2(W, #’) over all fermion 
self-energy parts, and Y*(f*) is the sum over all proper fermion self-energy 
parts. J7(f#) and J7*(f) are similarly defined. Square parts which are 
not quasi-proper are clearly redundant. Thus K(#, #, t°) is defined as 
the sum over all quasi-proper square parts and K*(#", ¢, #?) as the suni 
over all proper square parts. 
P2 
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To obtain the matrix element for a given process, one must consider 
the totally irreducible graphs for the process. ‘The corresponding integrals 
will be made up of factors Sp, 4p, y; and K(8,) (from the simple square 
parts). To make allowance for the insertion of self-energy parts into the 
lines of such graphs, the factors Sp and 4, must be replaced by Sp and 
A, which are implicitly defined by the equations (II. Eqs. (63), (64)), 


S4(2)=S,(H+8p(AZ*(A)Se(#),  . . . - (81) 
Ay (2) =Ap(*)+4—(A)*(2)45(8). . -.- . (82) 


These equations express the fact that a self-energy part is either proper 
or is a proper self-energy part attached by a single line to a self-energy 
part, which may be proper or improper. The insertion of all quasi-proper 
square parts in place of a simple square part is allowed for by replacing 
the factor K(S,) by K. Now a quasi-proper square part is either a proper 
square part or is a proper square part joined to a quasi-proper square part : 
by two meson lines, which may themselves contain self-energy parts. 
Thus, corresponding to (31) and (32), K is implicitly defined by the 
equation 
K(#, 2, A)=K*(A, 2, @) 

+fK*(#, #, t)4,(H)4,(A+P—)K(t, A+ —t, 8) d4t.  . (33) 


No infinity will arise from the integral over ¢t. A(#, #) is‘summed over 
all proper vertex parts only, and the effect of inserting vertex parts. 
into the totally irreducible graphs is allowed for by replacing y; by 
I'(@#, t?), which is defined by 
LEC, A=ys TAG). Sie oa ee 

All the infinities in the matrix elements will appear in the functions. 
Sp, 4;, and K. It must now be shown how these can be removed. 

Consider for the moment only totally irreducible parts. If S is a totally 
irreducible square part, then K(S, #, #, ¢?) contains divergent factors K(S,) 
from the simple square parts contained in 8. These simple square parts. 
can each be replaced by a single vertex at which 6A is operating, and the 
first term in the expansion of 6A (of order f*) can be chosen to cancel 
these infinities. (The higher-order terms in 5A will be defined later to. 
cancel certain other infinities which arise when irreducible square parts 
are considered.) With this understanding, the integral K(S, #, 2, #), as 
a whole, still diverges logarithmically. By splitting off the part obtained 
by setting 1 —f—t?—0 and summing over all totally irreducible square. 
parts, 


K*(#, 2, 8)=A"+i8A+K*(#, 2,8), . . . . (35): 


where A’’ is an infinite constant and K* is a uniquely defined finite 
function. With the same understanding about simple square parts, 
the logarithmically divergent A(V, ¢, @) can be split into a finite constant 
and a finite part. By summing over all totally irreducible vertex parts, V, 
we have 


A(t, 8) s Leg el F y ap eee-aieae 
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Totally. irreducible self-energy parts do not contain any square parts 
(see below), and can be split up as in electrodynamics (II. Eqs. (67, (77)) 


3 


2* (0) =A — 273K +BUtiy,,—iKo) + (Hy, —tko)S-(#), . . (37) 
IT*(t!) = A’ idx? +C((2)? +47) + ((A)2+e2)4,(2). 2. (38) 


L, A, B, A’ and C are infinite constants and 4,(¢), 8,(f) and A(#@, #) are 
defined by absolutely convergent integrals. Substituting these expressions 
into (31), (32), (33) and (34), equations are obtained for Si(t#), A/(t), 
K(@, ?, #) and I(@, t). S},, 44, [, and K, are defined to be the functions 
which are obtained if the infinite constants are simply dropped from these 
equations. Note that so far no account has been taken of self-energy, 
vertex and square parts which are not totally irreducible. 

It will now be shown that, by successive approximation, it is possible 
to develop power series in the coupling constant for finite functions Sj, 
Ay, K, and I’, even when graphs which are not totally irreducible are 
considered. Suppose that the series for these functions have been 
developed up to /*”. To develop the next term of the series for K*, 
consider the contribution from a totally irreducible Square graph, S, 
of order 1, (J<2n+4). The integral corresponding to this graph will 
contain factors S;, 4, and y; and factors K(S,) for the simple square parts. 
The exact integral is obtained by replacing these by 8,,, 4;,, (1, and K, 
respectively. These will not introduce any divergencies and will determine 
the integrand up to the order f2"+", The whole integral will itself diverge 
logarithmically, but the finite part can be separated from an infinite 
constant by the usual procedure to determine K,(S) and A’’(S) to the order 
St. Since 1>4, by summing over all such graphs, 8, power series in f 
can be developed for K, and A” up to the order f*”*4. 

Toobtain I’, to order f2"+?,sum inthe same way over all totallyirreducible 
vertex graphs, V, after substituting the corrected functions for the 
line, vertex and simple square factors. The logarithmically divergent 
integrals can again be split into an infinite constant and a convergent 
part to give power series for L and A, up to f*"**. Hence I’ can be found 
to the same order and I’, by dropping the infinite constant L. 

To obtain similar expressions for S},, and Ay,, it is not necessary to 
consider explicitly the inclusion of square parts; because the only 
irreducible self-energy parts are loops with two vertices, (II. Figure 5) ; 
and square parts can only be included as sub-parts of included vertex 
parts. They will thus be allowed for by the substitution of 1’, for ys. 
The development of S;, and 4}, will thus be just the same as in 
‘electrodynamics. 

Since the infinite constant A’ can be developed as a power series in f, 
the constant 5A can be chosen to satisfy 


Geo Oe ae A gene. (39) 


{This is consistent with the choice of the first term of 6A made above.) 
Then no infinities will arise from square parts provided the sum is taken 
over all such parts at each stage: K* is equal to K% and K, is identical 
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with K. But such summations always are taken in the procedure laid down. 
by Dyson. Thus, if this choice of 6A is made, no infinities will arise from. 
square parts if Dyson’s procedure for electrodynamics is followed for this 
interaction. In particular Dyson’s proof (IT. § 7), that the dropping of the- 
other infinite constants is equivalent to a ‘ charge > renormalization, 
if S=Zy GLE fe i Oe ee ee ae 
can be taken over directly, with the trivial addition that 5x? must be 
chosen to satisfy 
SOS (ela tA ee ee 
The renormalized ‘ charge’ f, can be equated to the finite experimental 
value and a finite S-matrix is thus obtained for effects arising from this 
interaction. Strictly speaking the elements of the S-matrix have been 
derived as a power series in f with finite terms. The convergence of the 
series itself has not been discussed. 
Exactly the same argument applies to the scalar interaction of scalar 
mesons with nucleons if y; in the above is replaced by —176,,, where 
6,, is a Kronecker 6-symbol operating on the spinor suffixes of and #. 


lof 


$5. OTHER MESON-NUCLEON INTERACTIONS. 

In all the remaining meson-nucleon interactions the condition of § 3 is 
not satisfied, and the factor corresponding to a meson line in an integral M 
will contain a double derivative of a 4,,-function. Thus, for considerations 
of convergence, a meson line will contribute to M a factor 


(Za;pi,)(Sa,p!)| ((Za,p')?+x2). 
The new convergence function K” is defined 
K?=F,--4(F—n- 1). ye ee 
This function has the same properties with regard to the convergence 
of the integrals now under consideration as K had for the interactions 
previously considered. The expressions (23) and (24) for F; and F are 
still valid so that 


K? — (3/2)K,+2E,—n+m,+n,—4.° . . . . (48), 
A primitive graph certainly converges if 
KP bey cute! tdci rie cabaneas 
but primitive graphs will be primitive derivatives if, by (43), 
(3/2) Riese ine Werte, in trel.  cemeet rs ere 


Owing to the appearance of » on the right-hand side of this inequality, 
there is no restriction on the number and degree of the primitive: 
divergents. An irreducible divergent graph may be found for any real 
process provided sufficient virtual processes are included. The integrals 
will only converge if the coefficients of all the divergent powers of the. 
basic energy-momenta vanish identically. That this is not always so: 
can be proved by any particular example. Thus for these interactions 
the renormalization of mass and charge is quite inadequate to remove: 
all the divergences. 
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The only exception to this is in the vector interaction of neutral vector 
mesons, in which the derivatives in the factor (17) come from the 
commutation relations themselves (16) and not from derivatives of the 
commutation relations as in (15). If the Stiickelberg (1938) form is 
used, in which two subsidiary fields are introduced, it can be seen that 
the derivatives arise from a vector interaction of the subsidiary scalar 
field. This can be eliminated by a transformation, and Dyson’s treatment 
of electrodynamics can then be applied without the necessity of intro- 
ducing a 6A term. This interaction has already been dealt with elsewhere 
(Matthews 1949 b), where it is shown that (17) can be replaced by 


{P(¢.,(z), $,(/))) >= 3%08,,4p(@—y).— . . . .* (46) 


§ 6. COMPARISON WITH RADIATION DAMPING. 

It is interesting to see how the renormalizing procedure when applied 
to mesons compares with the much more drastic subtraction procedure 
of Heitler’s (1941) theory of radiation damping *. Heitler’s theory has 
been most successful when applied to the scattering and creation of mesons 
by nucleons. In the cases where ordinary perturbation theory gives a 
cross-section which increases with primary energy, the inclusion of Heitler’s 
damping term gives a cross-section which decreases with primary energy. 
However, this situation does not arise in the scalar or speudo-scalar 
meson scalar interactions, nor for the vector interaction of neutral vector 
mesons. (If it did occur, it would be only for longitudinal mesons, and 
Feynman (1949 b) has proved that the cross-section is zero in this case. 
The proof is given for photons and is applicable to neutral vector mesons. 
However, it depends on the cancellation of certain terms, some of which 
do not occur if the mesons are charged. But for charged vector mesons 
Dyson’s procedure cannot be effectively applied.) Thus the renormalizing 
procedure is inadequate in all cases where an interesting comparison 
might have been made, and Dyson’s theory, as here applied to meson- 
nucleon interaction, is not, as might have been hoped, in any important 
case an effective and preferable alternative to Heitler’s. 
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XVII. CORRESPONDENCE 


Observation of Unimolecular Growth Steps on Crystal Surfaces. 


By L. J. GRIFFIN, 


Department of Physics, Royal Holloway College *. 
| Received December 16, 1949. | 
[Plates VII & VIII.] 


Aw extensive line structure on prism faces of beryl crystals has been 
revealed by microscopic observation. Evidence is given below to show 
that these lines represent steps only one unit cell in height; if this is 
the case the height of the step should be 9-17 Ax cos 30=7-9 A). 

The appearance and behaviour of these lines are the same whether the 
face is examined in its untouched, natural state, or after the deposition 
of a highly reflecting thin film of silver on the face. The sole difference 
between the two cases is that the visibility of the lines is much higher 
when the surface has been silvered. These lines must therefore corres- 
pond to genuine geometric surface details. 

As nearly as one can tell, these lines follow contour lines of the surface 
topography. The lines are widely spaced where the surface is relatively 
flat, and crowd together where there is a steep gradient. The only evident 
departure from contouring is that occasionally a single line terminates 
at a point in the face, instead of forming a closed loop. These lines must 
therefore represent steps on the surface. If the steps are taken to be one 
unit cell high, the point of termination must be interpreted as the end 
of a screw dislocation. 

An exact survey of the surface topography by interferometric methods 
should permit of a direct measurement of the step height at each line. 
So far attempts to apply this method have only sufficed to show that the 
step height is less than 34 A. From the fact that no visible kink appears 
where an interference fringe crosses one of the lines, one may also infer 
that the step height is less than 20 A. 

Further evidence will be presented below, in the discussion of fig. 3, 
that the lines represent steps all of the same height. The evidence thus 
points towards these steps being one unit cell high. 

The observation under the microscope of such small steps was 
surprising, and so contrary to normal practice or prediction, that it was 
felt that further evidence was needed before the results could be reported. 
They were not, therefore, made public outside this Department. During 


* Communicated by the Author. ' 
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‘the Bristol University Summer School on the Theory of Crystal Growth 
in September 1949, the theory of the behaviour of the steps originating 
from screw dislocations was developed much more fully by Frank, and 
it was seen that the correspondence was so exact between the predictions 
of this theory, developed for unimolecular steps, and the observations 
on beryl crystals, that one was led to the unavoidable conclusion that the 
observed steps were, in fact, only one unit celi high. 

The technique used in observing these steps consisted in first coating 
the face to be examined with a film of silver of reflectivity some 70-90 
per cent, the silver being deposited under conditions normally used for 
multiple-beam interferometry (thermal evaporation in a vacuum of 
pressure about 107° mm.Hg.). The surface was then examined in 
reflexion using a normal type of metallurgical microscope, bright field 
illumination being used with the field iris stopped down to give a fairly 
narrow illuminating beam, the specimen being racked slightly off focus. 

By going slightly off focus when observing the surface, a phase dif- 
ference is introduced between the diffracted wave fronts, and a certain 
intensity distribution is produced defining the position of the step edges. 
The precise nature of the intensity distribution is, of course, dependent 
on the position of the focal plane of the objective, non-complementary 
intensity distributions being produced on opposite sides of the true focal 
plane. In general, it has been found best to work with the edge imaged 
by an intensity distribution consisting of a broad minimum between two 
narrow maxima (see figures). These unimolecular step edges have been 
observed with objectives of N.A. ranging between 0-058 and 1-3, the 
visibility of the effect apparently decreasing with decreasing numerical 
aperture of the objective ; it is, however, difficult to determine whether 
this decrease in visibility is primarily due to the decrease in the observing 
numerical aperture, considered as the collector of a cone of light of certain 
-angle, or is caused, partly or entirely, by other factors dependent on the 
value of the numerical aperture. 

The deposition of a reflecting silver layer, some 250-400 A.U. thick, 
greatly increases the visibility of the diffraction effect produced by the 
step edges. It does not seem to alter the nature of the intensity 
distribution in the image of the edges as obtained from an unsilvered face. 
"This was to be expected from the evidence obtained from multiple-beam 
interferometry (Tolansky 1948) that such silver films, when deposited 
under correct conditions, contour the surface of the substrate to within 
very close limits. 

It would seem probable that improvements could be made in the 
technique of observation which would further increase the sensitivity. 

The topography of prism faces of beryl crystals is in general complex, 
with a tendency to form occasional low plateaux of a typical “ razor- 
blade ”’ shape, similar in outline to the outer loop of fig. 2. On the surface 
-of two such main plateaux the two very simple cases shown in figs. ] and 2 
were found. 
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Fig. | (magnification 625 x ) shows the simplest case of a single screw 
dislocation, imaged here as a dot, from which runs a spiral edge, as. 
predicted by the theory (Frank 1949 a). The c-axis of the crystal is 
here, as in all the figures, parallel to the almost straight edges. The 
light area is a region of very coarse structure, the general, irregular, and 
randomly distributed, surface markings being probably due to the 
deposition of foreign material after the beryl crystal had ceased to grow. 
This irregular and random marking is of general occurrence on faces of 
beryl, varying considerably in density. 

Fig. 2 (magnification 1050 x) shows the next most simple case, that 
of growth from the step joining & right and a left handed screw dislocation 
(Frank 1949 b). According to the theory unclosed loops like the inner 
one, attached to a pair of dislocations, extend during growth of the 
crystal in such a way as to generate successive closed loops, exemplified 
by the outer one. (The triangular area cutting into the closed loop 
probably consists of an over growth of foreign crystalline material). 
The two dislocations in this case are also imaged as dots, this being 
the usual type of imaging on these prism faces. 

Figs. 3 and 4 show the more normal type of behaviour seen on beryl, 
the above simple cases being exceptional. In fig. 3 (magnification 625 x ), 
the steps advancing from the right are split into three sections by a 
succession of three crystal inclusions, which act as obstacles. Of the 
three sections, the central and lower sections first combine again to. 
form a single section. Counting from the last edge before division 
commences on the right, there are 37 edges in both the central and lower 
sections before junction occurs again. Similarly, considering the junction 
beyond the third obstacle, there are 57 edges in both the upper and 
lower sections between the first division of the edges and the first complete 
edge. It must, therefore, be that, in general, each step is a like unit 
step. For a step one unit cell high this type of behaviour is immediateiy 
understandable, and in fact necessary, but for steps of two or more unit 
cells it would appear very probable that in these circumstances the 
steps would either split into smaller units or combine to form multiple: 
steps. 

The slow closing in of the edges beyond any inclusion, and the rapid 
elimination of the re-entrant angle formed at the junctions of the 
clifferent sections, is typical. 

Fig. 4 (magnification 500 x ) shows the effect produced by a distribution 
of disclocations “ dominated’ by the steps originating from another 
group, far off the picture to the left. Some of its features are discussed 
by Frank, below. 

Similar behaviour to that reported above has been seen on two other 
crystals of beryl, in these cases of typical emerald appearance and habit, 
while the above crystal was of aquamarine type. The three crystals 
came from widely different localities (e.g. the Urals and California), so. 
that this would seem to be a general type of behaviour of beryl crystals... 
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Figs. 1, 3 and 4 negative. Fig. 2 positive. 
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Taking into account the shortcomings of two-dimensional nucleation 
theory (Burton, Cabrera and Frank 1949, Cabrera and Burton 1949), 
the general validity of the dislocation theory of crystal growth would 
seem to be established. Further confirmation of the theory should be 
obtainable not only by observations of similar static phenomena on other 
types of crystals, but also by observing successive small amounts of 
zrowth or solution on faces of crystals, the topography of which had 
been previously determined. By this means a kind of cinematograph 
picture could be built up of the behaviour of the steps. In the case of 
beryl, further growth is not possible, but it is hoped to observe solution 
effects and thereby test the theory further. 

Some small rod-like structures can be seen in the lower part of fig. 4. 
This rod type of structure is of common occurrence on the beryl crystals 
examined and occurs invariably on steep gradients running in the general 
direction of the c-axis. The rods are often found with one or more 
dislocations at their base and usually show a somewhat irregular step 
structure. 

The difference in density of occurrence of screw dislocations per unit 
area can be considerable, both between different areas of the same face, 
and erystallographically equivalent faces of the same crystal. Thus 
there are some 100 screw dislocations in an area of 7-6 X 1073 sq. mm.. 
of one face, on other parts of which areas of several square millimetres 
show no sign of dislocations. An equivalent prism face of the same 
crystal shows no sign of dislocations over an area of several square 
centimetres. There is a definite tendency for dislocations to be closely 
spaced along lines either roughly parallel or perpendicular to the c-axis, 
as one would expect from general dislocation theory. 
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Note on the Structure of a Crystal Surface. 


By F: C. Frank, 
H. H. Wills Physical Laboratory, Bristol *. 


| Received December 23, 1949. ] 


THE observations reported by Griffin (1950) provide an unexpected 
direct confirmation of the theoretical conclusion that dislocations play 
an essential role in the growth of crystals. Indeed, they provide the 
first direct visual evidence of the existence of dislocations themselves. 
They also offer promise of a much closer insight into the actual conditions 
under which the growth of natural crystals occurred. If it also turns 
out that a similar technique will reveal details of ** molecular ” thickness 
on other crystals, growing in the laboratory, the possibilities for enlarging 
our knowledge about crystal growth processes are, of course, enormous. 

The quantitative theory of crystal growth which has been developed 
in various papers by Burton, Cabrera and the writer (1949, 1950) (some 
of which, as yet unpublished, were presented at the Bristol Summer 
School on the Theory of Crystal Growth) has so far dealt mainly with 
crystals formed from spherical or cubic molecules, with growth occurring 
on crystal faces of high symmetry. In qualitative principle the theory 
applies equally, though with complications, to the (1010) face of Beryl, 
studied by Griffin. There appears to be no reason to doubt that the 
essential building unit of Beryl is the Sig 0, ring with its attendant two 
aluminium and three beryllium ions, making half the unit cell found by 
Bragg and West (1926). These units may be regarded as * pill-boxes ”, 
7-9 A. in diameter and 4-6 A. in height, close-packed side by side in the 
basal planes of the crystal, and one on top of the other in the direction 
of the hexagonal axis. In this direction alternate Si, 0,, rings differ in 
orientation by 25°. This alternation secures the best coordination of 
oxygens about the cations, and doubles the X-ray unit cell, but may be 
considered a secondary detail of the building principles so far as we are 
concerned. On a (1010) face these “ pill-boxes” lie on their sides 
packed in flat layers in rectangular array, and growth would be expected 
to proceed by the spreading of such ‘‘ monomolecular’? layers by 
accretion of new units at their edges. At the low supersaturations 
essential for the growth of good crystals such layers would not develop 
on a free flat surface except from the points of emergence of dislocations 
with a component of the displacement vector normal to the surface 
(screw dislocations ” for short, though it is not essential that they 
should be pure screw dislocations). The way in which the edges of 
monolayers develop as ‘ growth spirals’ centred on the dislocations 
have been described elsewhere (Frank 1949), and beautifully illustrated 
in Griffin’s figs. | and 2. 


*Communicated by the Author. 
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~The dependence of the two-dimensional crystal growth rate on:the 
orientation of the step line imposes a characteristic distortion of the. 
growth spirals and loops, exhibiting the symmetry of the crystal face. 
A fairly simple theory can be given to account for the shape assumed 
by each growth layer as it expands, the characteristic ‘“razor-blade ”’ 
Shape seen on the beryl crystal being a typical example. In essence 
it is the same theory as was developed earlier (see for example, Friedel 
(1926)) for three dimensional crystal growth on the assumption that the 
growth rate at any point on a crystal surface depends only on’ the local 
crystallographic orientation of the surface, and its environment. We 
now appreciate that this is a far from correct assumption in that case, 
but in two dimensions it should be much more nearly true. In the 
latter case, a straight edge corresponds to a sharp minimum in the growth 
rate as a function of orientation. Truly sharp minima, and therefore 
truly straight boundaries are not in general to be expected. When the 
rate of advance, v, of the step, depends on its orientation % according 
to the law 


v=, (1+ By), 


we can show that the resulting boundary will have a radius of curvature 
(1 + 2B) times its distance from the centre of the figure. A simple 
theory, applicable in certain circumstances which are particularly likely 
to occur in growth from solution, makes the rate of advance v proportional 
to the number of © Frenkel kinks ” in the step line. Then it follows that 
B=e”"T/8, where ¢ is a neighbour-neighbour interaction energy 
which we may take to be } of the latent heat of solution. In this case 
the nearly straight edge of the razor-blade shape in Griffin’s figure (2) 
has an apparent radius of curvature 100 times its distance from the 
centre, though it may be more (?.e. the apparent curvature is possibly 
due to optical distortion). Thus ¢/kT is 6, or more. By contrast, 
the growth figures on a synthetic emerald grown hydrothermally by 
Nacken (1948) (Photographs by Griffin, not reproduced here) are 
similar, but have strongly curved edges to the “ razor-blade ” indicating 
B=8, whence $¢/kT=4-16. Supposing the natural and synthetic 
erystals grew from similar media, we may deduce that the absolute 
temperature of growth of the natural crystal was only about two-thirds 
Nacken’s working temperature. Information of this is unfortunately 
lacking, but the indication that the natural crystal probably grew at 
a moderate temperature may be of interest. 

The two ends of the razor-blade are very nearly ares of the same 
circle both for. natural beryl and Nacken’s crystals, indicating that the 
growth rate is practically independent of orientation for the range of 
orientations occurring in them (-+-30° from the normal to the hexagonal 
axis in the case of the natural crystal). 

Fig. 5* shows polar diagrams of two dimensional growth rate which. 


~ 


* Figures in this letter are numbered from 5 onwards to simplify reference 
to figures in Griffin’s article, figs. 1 to 4. 
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would account for the observed two-dimensional growth figures. These 
polar diagrams are naturally uncertain for the range of orientations 
which do not appear in the figures: in such ranges one can only say for 
certain that the polar diagram lies outside the circles whose diameters 
_are the lines connecting the centre to the corners of the figure. 


Fig. 5. 
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The “razor-blade” outline of a growth layer on the (1010) face of beryl 
(hexagonal axis horizontal) and the associated polar diagram of rate of 
erowth as a function of boundary orientation, (a2) for the natural crystal, 


(b) for Nacken’s synthetic crystal. The form of the polar diagram is 

. . =) : 
uncertain where represented by a dashed line and coincides with the 
‘“* razor-blade ”’ outline on circular ares at the ends. 


Observation of the spacing between turns of a growth spiral, or of its 
form where it joins the dislocation, enables us to measure a characteristic 
length, “the size of the critical nucleus’, whose order of magnitude is 
ap/kTinwn where a is a lattice spacing of the crystal and « is 
the saturation ratio. From fig. 2 this length is found to be about 1 micron. 
whence ¢/kTin« is about 10°. In conjunction with the former estimate 
of ¢/kT we may deduce that the degree of supersaturation at which 
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‘growth took place was about | per cent. These estimates are not intended 
to be final or accurate, but are meant to indicate the kind of information 
which we should be able to derive from the type of observation Griffin 
has been able to make. . 

Whereas Griffin’s figs. (1) and (2) show the simplest kind of growth 
initiation centres, he has other photographs showing initiation by larger 
groups of dislocations, such as occur in his fig. 4, though here they are 
not acting as an initiating group because they are “dominated” by a 


Fig. 6. 
\ 


Three successive stages in the passage of a dominating sequence of growth 
steps over a dislocation. Each line represents a unimolecular step, down 
to the right. The next stage in this succession will resemble the first. 


Schematic perspective sketch showing crystal surface with two unimolecular 

growth steps and a short shallow slip zone. The terminations of the 
latter provide three left-handed and three right-handed dislocations. 
Further growth (dominated by the step sequence) wil! produce a pattern 
corresponding to Griffin’s fig. 4. 


more active group lying off the picture to the left, in the way described 
by Frank (1949, b). The sketches in fig. 6 show three successive stages 
in the passage of a dominating sequence of growth fronts across a single 
dislocation, stages which may be recognized at various places in fig. 4. 
Fig. 4 also shows, as well as the randomly distributed dislocations, two 
groups lying approximately on a straight line, making a closely spaced row 
of right-handed dislocations at the top of the picture, and of left handed 
dislocations near the bottom. This is the kind of distribution of screw 
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dislocations which results from the formation of a limited slip-zone, as. 
illustrated in fig. 7. One may readily see how subsequent growth on a 
surface of this configuration, under the domination of a close sequence 
of step-lines originating elsewhere, produces the general pattern of fig. 4. 

It is interesting to consider the origin of such slip zones, and similar 
arrangements of dislocations to be found elsewhere in pictures taken by. 
Griffin, in which, however, the line is parallel to the hexagonal axis, instead 
of transverse to it. A likely cause is the stress resulting from inhomogeneity 
of ionic replacement in the crystal. The fact that the crystal is coloured 
is sufficient evidence that ionic replacement occurs. A change in the 
composition of the solution will produce inhomogeneity in the crystal, 
which must evidently lead to lattice curvature. Now, consideration 
of the growth of a curved crystal from a mixed solution shows that it 
will suffer oscillations in the lattice spacing, with a wave length of the 


Fig. 8. 


(2) 


(a) A form of growth front occasionally observed ; (6) the normal form of growth 
spiral (i) is transformed into the shape (iii), with reversed curvature, as 
in 8(a), by a slight dissolution following normal growth. 


order of magnitude of its shortest dimension, accompanied by considerable 
stresses. Dislocations already present may move under these stresses, 
but, additionally, the surface stress in new growing layers (particularly 
shear stress and tension) can, relatively easily, produce surface faults, 
which become true dislocations when buried under additional layers. 
This is probably the chief source of new dislocations, lineage boundaries 
and the like, in slowly growing crystals. Of course, the initial dislocations 
needed to make the first crystal seed capable of growth are generated with 
ease at the high supersaturation needed for nucleation. 

On certain crystals, Griffin finds that the growth fronts have a reversed 
curvature near their terminations on dislocations, as in fig. 8(a). This 
appears to indicate that the last phase in the development of the crystal 
was a slight dissolution, before all change ceased. Fig. 8(b) shows how the 


‘ 
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normal growth spiral (i), retracting in dissolution by way of stage (ii), 
proceeds to wind the spiral the other way, close to the dislocation, as in 
stage (iii) which corresponds to 8 (a). 
H. H. Wills Physical Laboratory, F, C. FRANK. 
_ The University, Bristol. 
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Transfer Rates for Surface Films of Hell. 


By E. J. Buree and L. C. Jackson, 
H. H. Wills Physical Laboratory, University of Bristol *. 


PUBLISHED measurements of the rate of transfer of Hell by the surface 
film fall into two groups. On the one hand there are the pioneer studies 
of Daunt and Mendelssohn (1939) which showed that the rate depends 
only on the smallest periphery above the liquid helium level and on the 
temperature, increasing from zero at the A-point to a flat maximum 
of about 7:-5<107> cm#/sec/em. at 15° K. These observations were 
substantially confirmed by the work of Webber, Fairbank and Lane 
(1949) and the later work of Atkins (1949). 

On the other hand Atkins in his earlier work (Atkins 1948 a) observed 
rates up to ten times those of Daunt and Mendelssohn and showing a 
marked dependence on the distance of the higher helium level from the 
rim of the beakers but independent of the distance between the two 
liquid levels. He used a capillary beaker terminating in a spherical bulb 
and the results were attributed to temperature gradients and possibly 
the geometry of the vessel. 

Similarly de Haas and van den Berg (1949) have recorded high creep 
rates and a dependence of the rate on the difference between the helium 
levels. They found the rates of emptying the beaker were those to be 
expected for a syphoning process, 7.¢. dependent on the difference of 
height of the helium levels. This type of process will not, however, explain 
the rates of filling initially observed by Atkins, and confirmed later at 
Oxford and at Bristol (for contaminated surfaces), which were independent 
of the difference of height of the levels but dependent on the height of 
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the rim of the beaker above the highest helium level. Their experiments 
were carried out in the absence of radiation and they attribute the 
difference between their results and those of Daunt and Mendelssohn to 
this factor. 

Bowers: and Mendelssohn (1949) showed that absence of radiation had - 
no effect on the normal creep rates, and that the high creep rates observed 
by Atkins in his first experiments were due to films of solidified gases 
(c.g. air) on the walls of the vessels, In his later experiments Atkins has 
confirmed the correctness of this explanation. It may be of interest to 
record our own observations on creep rates collected during the course of 
our experiments on helium II films (Jackson & Burge, 1949) and designed 
to test the similarity of our experimental conditions to those of Daunt 
and Mendelssohn in their measurements of creep rates and film thickness. 
These observations were made with a cylindrical beaker enclosed in the 
combined uniform temperature enclosure and radiation shield shown in 


Fig. 1 (a) Fig. I (0) 


fig. | (a). The walls are silvered except for two narrow longitudinal slits. 
A half-watt neon lamp with a filter of Chance’s ON19 glass served as the 
source of illumination. The creep rates agreed closely with those of 
Daunt and Mendelssohn. 

After the publication of Atkins’ first note (Atkins, 1948 ay the cylindrical 
beaker of bee 1 (a) was replaced by one consisting of a glass capillary tube, 
length 2-2cm., internal diameter, 0-04 em. terminating in a nent 
bulb of aietniseee 0-8 cm. fig. 1 (b). The creep rates observed again agreed 
very well with those of Daunt and Mendelssohn. The appearance of the 
note of Bowers and Mendelssohn led us to repeat our observations 
(1) with the ‘* clean ”’ Atkins type beaker, and (2) after admitting some 
50 cm? of air to the cryostat, a small portion of which condensed upon the 
beaker. The creep rate with the “ clean’ beaker was again normal but 
the maximum rate with the contaminated beaker was some seven times 
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greater and a marked variation with height was observed (tig taede: 
T=1-90° K). We confirmed Atkins’ observations that the filling Toten 
depended on the height of the rim above the highest helium ia and 
not on the difference of the levels. In some of the experiments with 
the contaminated beaker the observations were continued until the 
levels of the liquid helium inside and outside the beaker were the same. 
The inner level was then found to overshoot and the oscillations first 
described by Allen and Misener (1939) and studied in detail by Atkins 
(1948 b) were observed. In one experiment twelve complete periods were 
timed, giving an average value for the period of oscillation of 10-4 secs. 


Fig. 2. 
Creep rate (cc/sec/cm x 10°) 


o Emptying. 


O Falling. 


Clean beaker. 
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Assuming the film thickness to be inversely proportional to the height 
and inserting the above value of the period of oscillation in the formula 
given by Atkins, the film thickness at a height of 1 cm. was calculated. 
The value obtained, d=11-16 X 107& cm. (T=1-07° K) agrees well with 
the early results of Atkins. Unfortunately similar observations had not 
been made with the clean beaker. 
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* 
XVIIL. Notices of New Books and Periodicals received, 


Introduction to Atomic Physics. By 8. Totansky. Third Edition. [Pp. 371.] 

(Longmans Green.) Price 15s. 

THE appearance of this text for the degree student is welcome because the 
second edition with its 1946 appendix was already somewhat out of date. 

New sections on recent developments have been included in various subjects. 
The chapters on cosmic rays and the cyclotron type of machine have been 
brought more up to date and that on the nucleus has been enlarged. Of the 
two new chapters added, one deals entirely with the neutron and its properties ; 
the second concerns nuclear fission and atomic energy. 

It is in one way unfortunate for the author that some of the subjects dealt 
with are developing at so remarkable a rate that the information is already 
out of date, e.g., in the sections dealing with mesons and the photographic study 
of nuclear and electron tracks. ’ 

In our experience students have been almost unanimous in praising the second 
edition as a stimulating source of valuable information. They will continue to 
do so with the new edition. The work is relatively free from the faults that 
marred the first edition, though the reviewer notices one unfortunate statement 
on page 318 in the paragraph on the betatron, namely, “The mass increase 
limits the ultimate velocity ’’. It is feared that the student will take this to 
mean that relativistic change of mass sets a limit to the energy obtainable, 
which, of course, is untrue. i H. H. P. 


Ions, Electrons and Ionizing Radiations. By J. A. CRowTHEeR. (8th Edition.) 
[Pp. 322.] (Edward Arnold.) Price 21s. Od. 


THat this work has been accepted for many years as a standard text-book 
is shown by the fact that the demand for it has required an eighth edition. 
One reason is that it is eminently readable and contains no information which 
is superficial or unreliable. The main changes which have been introduced 
since the last edition, eleven years ago, have naturally been in the field of 
nuclear physics and atomic energy. In this field there is such a wide choice of 
material that there is room for individual preference. The reviewer would like 
to have seen some reference to Powell’s development of the photographic 
method for the quantitative study of nuclear particles. It is now ten years. 
since he pointed out the potentialities of this method, which are now being so 
efiectively exploited as to speed up the tempo of advance many times. On the 
other hand, one can welcome the author’s refusal to throw on the scrap-heap 
descriptions of early discoveries in the vain hope of circumscribing completely 
the subject of atomic physics within the limits of a single volume. This point 
of view will be shared by many teachers, particularly those who, like the 
reviewer, have watched the subject develop through the past fifty years. 


A.M. T. 


Faraday’s Discovery of Electromagnetic Induction. By T. Marrty. [Pp. 160. 
(London : Arnold.) Price 9s. 


THe author of this very readable account of a great discovery and the events 
which led up to it is the editor of Faraday’s diary published before the war. 
In a connected story, illustrated by excerpts and drawings from the diary 
readers can catch some of the atmosphere of scientific investigation a century 


and more ago and sense the working of mind and intuition of a great 
experimenter, 


[The Kditors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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